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ABSTRACT 


A generalized computer programme has been deve- 
loped to optimize the insulation thickness, sink and 
dimensions of cold storage. The total cost comprising 
the initial and running costs based on the present worth 
method of economic analysis has been taken as the objec- 
tive. Optimum insulation thickness has been determined 
for the ceiling and that part of the walls which is 
above and below ground level. The optimum sink to which 
the structure must be built, in order to minimize func- 
tional energy and cost, has also been determined. 

Optimum length, breadth and height of the room 
have been determined for a given volume of the storage 
space. 

Powell's method of unconstrained minimization 
coupled with the interior penalty function method to 
convert the constrained problem to the unconstrained one 
and the quadratic interpolation method are used in the 
optimization procedure. 

Optimum insulation thickness has been deter- 
mined for two types of insulating materials; Thermocole 
and Polystyrene foam insulation. 
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Optimum values of the design variables and the 
cost of the cold storage have been appraised for both 
cooling and freezing. 

The vapour-compression refrigeration system 
using Ammonia as the working fluid has been selected in 
computing the cost of power needed for refrigeration. 

Threlkeld # s classical approach involving the 
concept of sol-air temperature has been made use of in 
calculating the cooling load on the system. 

The variation of total cost with insulation 
thickness and sink has been studied for different dimen- 
sions/ inside temperatures and different insulating 
materials. 

Optimum values of the parameters determined by 
the hourly load calculation method has been compared 
with those obtained on the basis of standard design 
conditions . 

The effect of orientation of the building/ on 
optimum values of the parameters and the cost/has been 
studied. 

The decrease in functional energy and the 
increase in total cost is studied for suboptimal 
solutions. 
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I NTRODUCTION 

1 . 1 DESCRIPTION 

India, being mainly an agricultural country, 
produces large quantities of agricultural products of 
different varities. Some 20 to 30/ of these products 
are wasted in the areas where they are grown, because 
of unavailability of a controlled environment. The 
estimated loss is about Rs. 1500 crores per year due to 
the perishing of these products Ll]. 

Of the many perishable products the preservation 
of potatoes is the major application covering 30/ of the 
cold stores in India [ 2 ]. The use of cold stores has 
been further diversified to store other commodities like 
fish and dairy products (5/), fruits, vegetables, flower, 
jaggery and other items (15/). 

Cold storages maintain a low temperature and 
suitable humidity to prevent the spoilage of perishable 
products in one region and make them available in the 
off season as well as in other regions where they are 
not harvested. The growers also get a due share of the 



profit by way of not selling their seasonal produce at 
throw away prices. 

The rapid developments that are talcing place in 
the cold storage industry/ has led to greater importance 
being paid to the design aspect. Unlike the arbitrary 
design of earlier years, new vistas have been opened up 
in cold storage designs. The emphasis has been shifted 
to constructing cold storages and refrigerated warehouses 
with larger volumes, capable of storing a larger amount 
of commodities at a lesser cost. Improvements are being 
made in this field to minimise the functional energy 
requirements and the costs wherever possible. 

Generally, a number of alternatives can be 
thought of for designing a new system or improving an 
existing one. In cold storages the structural cooling 
load, which is dependent on environment, is quite subs- 
tantial apart from the commodity cooling load. The 
environmental parameters contributing to the structural 
cooling load are probabilistic and uncontrollable. Hence 
parameters like the thermal insulation, the area of the 
walls and the roof exposed to solar radiation and external 
environment, the orientation of the building, selection 
of design conditions and data, and the characteristics 
of the refrigerating system used for cooling are to be 
closely controlled in order to ensure optimum functional 
energy at lowest cost. 
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In this regard, energy saving can be achieved 

by decreasing the solar radiation falling on the walls 

and the roof of the cold store by sinking the structure 

partly or fully below the ground level, exposing lesser 

surface area to solar radiation. This also results in 

decreasing the inf lit. ration heat load owing to the foot 

that heavier air stays in that part of the structure 

which is below ground level. Another advantage of sirvk- 

is that the underground 

ing the structure partly or wholly under-groune^ tempera- 
ture of the soil is lower than the ambient temperature 
at any time of the day and moreover the underground 
temperature fluctuations with time is not as drastic as 
that of ambient temperature. Hence relatively more 
stable conditions prevail in partially underground storages 
imposing lesser transient load on the refrigeration system. 
Further, this system would serve a better purpose in 
places where electric breakdowns are frequent. 

It is also seen that the insulation requirement 
for underground: cold storages is considerably lower than 
that for cold storages which are above ground ievel, as 
a result of lesser heat transmission through underground 
walls. 

The concept of an underground cold storage 
however has some disadvantages; like extra digging and 
underground constructional cost and unforeseen problems 
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due to seepage of water from the underground. 

The shape and orientation of the structure also 
Plays a vital role in energy saving, since the west and 
east facing walls are subjected to longer exposures of 
solar radiation, the area of these walls should preferably 
be lesser than those of the north and south facing walls. 
Other orientations can be resorted to if it results in 
decreased structural cooling load and lesser costs. The 
best shape of the cold storage is that which renders 
minimum surface area for a given volume, rendering reduced 
cooling load and lesser insulation requirement. Ohe 
shape of the cold storage is influenced by other prac- 
tical problems like the constructional cost, maintenance 
and handling etc* 

Lastly, energy saving can also be achieved by 
selecting proper refrigerating machinery and regular 
maintenance in order to achieve a high coefficient of 
performance. 

In essence the optimum cold storage design is 
envisaged as having both financial feasibility and that 
functional energy which renders the overall cost for 
the cold storage a minimum. 
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1.2 LITERATURE REVIEWS 

The work available in the field of design of 
refrigeration and air-conditioning systems pertaining to 
cold storages and refrigerated warehouses is summarized 
in the following paragraphs. 

1.2.1 ANALYSIS OP COOLING LOAD CALCULATIONS AND 
HEAT TRANSFER THROUGH WALLS 

The analysis of cooling load involves parameters 
such as the ambient air temperature/ direct and diffused 
solar radiations/ air velocity/ characteristic of the 
enclosure walls and the orientation of the enclosure. 

Threlkeld L 3] has extensively described the 
process of heat transmission occurring through solid 
boundaries of a structure when a temperature differential 
exists between the internal and external environments. 

He has outlined a detailed procedure for the computation 
of solar irradiation on the walls and roof. Periodic 
heat transfer through the walls has been explicitly dealt 
with through the concept of sol-air temperature. His 
analysis essentially deals with heat transfer through 
homogeneous walls. 

Mackey and Wright [ 4 ] have described In 
detail the periodic heat transfer through composite 
walls or roofs. They have also proposed conversion 
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equations to reduce a composite wall problem into that 
of an equivalent homogeneous wall. The following three 
methods have been developed to account for the periodicity 
of external conditions: 

(i) Threlkeld's classical approach [3] 

(ii) Transfer function method [5] 

(iii) Finite difference method [s]. 

Threlkeld's classical approach involves the 
concept of sol— air temperature and is used in the pre- 
sent work. 

In the transfer function method the various 
components of space heat gain are added together to 
get an instantaneous total rate of space heat gain, 
which is then converted into an instantaneous space 
cooling load through the use of weighting factors called 
coefficients of 'room transfer functions'. The transfer 
function is nothing more than a set of coefficients that 
relate an output function at a given time to the value 
of one or more driving functions at the time and previous 
times. 

Kadambi and Hutchinson [ 6 ] have described an 
approximate technique to determine one dimensional 
transient heat transfer through walls and roofs, in the 
form of the Finite Difference method. The basic 



7 


simplicity of approximate method contrasted with analy- 
tical techniques is asserted in their work. 

Probably, the best and most relevant descrip- 
tion of cooling load calculations, from the point of 
view of practical designs, is outlined in ASHRAE, Hand- 
book of Fundamentals [ 7 ]. The methodology and equations 
for hour-by-hour load calculations, used in the present 
work, is outlined in the handbook. 

1.2.2 ECONOMICAL DESIGN AND OPTIMIZATION 

With the limited availability of energy and its 
consequent increased cost, it has been necessary to use 
energy optimally. Extensive work has been done in mini- 
mizing the functional energy and cost of cold stores and 
refrigerated warehouses which are above ground level. 

M cClur e 03 has described a method for the 
optimum design of building systems to reduce energy 
requirements. 

Bonar [9] has analysed the different parameters 
and factors that affect the economics of a refrigeration 
system for a warehouse with reference to the minimiza- 
tion of total costs. Thermal insulation is found to be 
the foremost factor that affects the functional energy 
requirement and the energy cost as well. 
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Prasad [ 2 ] has provided a break up of the total 
cost of a cold storage as follows: 

building work (30/), steelwork and cladding 
(13/), refrigeration plant and electric (21/) 
and insulation ( 31 /). Thus it can be seen 
that insulation alone accounts for about one 
thirds of the total cost. This gives a statis- 
tical picture as to the method of controlling 
the total overall cost of the whole system. 

Spielvogel [lo] has presented an interesting 
article as to how the use of more insulation than re- 
quired can increase the energy consumption and cost. 

Gupta [ll] has done extensive work on automated 
optimum design of refrigerated warehouses and air condi- 
tioned buildings. He has formulated the design problem 
as a nonlinear mathematical programming problem and has 
used multidimensional optimization techniques to solve 
it. His work presents a probabilistic optimum design 
procedure by considering the randomness of the input 
parameters. He has used the approximate partial deriva- 
tive method to evaluate the cooling or heating load. 

His. work is an example of the application of extremal 
distributions in the optimum design of refrigeration 
and air conditioning systems. 
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Heinze [12] 

has made a comparative study of the 
single and multistoreyed cold storages, and the depen- 
dance of insulation thickness on storeyed constructions. 
Single storey cold storages are suitable for places 
where the cost of land is quite cheap, whereas multi- 
storey cold storages are suitable to localities where 
land and labour are very costly. In such cases multi- 
storey structures has been found to be more economical 
provided mechanical handling is used. 

The work of Claesson and Holnquist [ ! 3 ] envisages 
the storage of cold and frozen food in underground mined 
rock caverns. Their work gives an idea of the energy 
saving from underground storage. They have suggested 
that such storages can be preferred where disused mines 
are readily available for underground construction. 

Prasad et al [ ll], in an analytical study, have 
determined the expressions for optimum insulation thick- 
ness and sink of an underground cold storage in terms 
of various non-dimensional parameters. He has determined 
a quadratic relationship between the approximate optimum 
sink and insulation thickness. This paper forms the 
basis for the present work which determines optimum 

parameters, for a cold storage of any shape through 

* . 

iterative numerical method. 
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1.2.3 COMPUTER APPLICATIONS 

Some years ago -a any extensive design analysis 
of air conditioning systems/ in which computations were 
performed manually, was considered impractical. But with 
the advent of computers in the field of design such 
extensive design analysis has become very common. 

ASHRAE [ 7 ] gives a methodology and equations 
for constructing computerized routines to duplicate 
results obta ined with the tabular data, and to examine 
the effect of various systems and operating schedules on 
the space cooling load. Even though the speed of norrrpti— 
tations performed by a modern digital computer enables 
a rigorous approach to be formulated for calculating 
cooling loads, there are very few computer programs in 
use today where exact cooling loads are calculated. 

Lokmanhekin and Henning er [id] have developed 
a comprehensive program that establishes cooling load 
requirements, energy estimating , equipment sizing, system 
simulation and economic analysis. Their program uses 
the advanced method of convolution principle to account 
for the thermal storage effect of a building. They 
have concluded that the convolution principle outlined 
in [ 16 ] gives more realistic values of peak loads and 
associated times of occurrence* An interesting feature 
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of their work is the thermal load plot which yields plots 
of the load profile of any space for any period of time. 
Comparison between plots permit grouping of compatible 
spaces into system control zones. 

Evers [l7] has developed a series of three com- 
puter programmes called the E cube series. These pro- 
grammes compute the energy requirements/ equipment selec- 
tion and energy consumption and economic comparison. The 
energy requirement programme takes design point values 
for seven components of thermal load and the base com- 
ponent electrical load and distributes them over each 
hour of the year in accordance with dry-bulb and dew 
point temperature variations. The programme further 
evaluates the effects of thermostat setback or periodic 
system shutdown and thermal storage of any magnitude 
including temoerature lag effects. 

The equipment selection and energy consumption 
program determines the actual energy consumed by the 
various pieces of equipment used to meet the hourly 
requirements. The economic comparison programme is 
oriented towards making decisions concerning investments/ 
annual operating costs and methods of financing. 



1.3 PRESENT STUDY 
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cost per year which is the objective function in the 
optimization problem, is calculated. The objective func- 
tion is minimized to achieve optimum values of the 
design parameters. The constraints on the design para- 
meters are formulated on the basis of practical considers' 
tions. Both one dimensional graphic and mathematical 
programming techniques have been used to obtain optimum 
values . 

In general the contributions of the present work 
can be stated as follows - 

1) Development of automated optimization procedure 
for the design of cold storages. 

2 ) Application of hourly weather data in the 
optimum design of such systems. 

3) Analysis of the effect of sinking the structure 
below ground level . 

4) Comparison of suitability of two different 
insulating materials. 

5) Comparison of optimum parameters for two 
different orientations of the building. 

6) Comparison of optimum parameters obtained by 
load calculations based on peak design condi- 
tions with those obtained by hourly load 
calculation. 
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In essence a generalizes computer programme 
has been developed for the opti™ design of cold storages 
and refrigerated warehouses. 

The cost of electricity for calculating power 
charges was based on data from standard sources. A 
logistic curve was fitted through these past data to 
forecast the costs for future, in order to determine 

core realistic power cost charges per year by present 
worth method. 


In order to generalise the program empirical 
relations for the properties of refrigerants, needed in 
the calculation of the coefficient of performance of the 
refrigeration system, are incorporated. These empirical 
relations depend only on the condenser and evaporator 
peratures or the system. The evaporator temperature 
selected on the requirements of inside temperature. 


Lastly the feasibility of using a suboptimal 
design was studied, having larger insulation of the 
order of 20 to 30/. 
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PROBLEM FORMULATION 

2.1 MATHEMATICAL MODELLING 

2.1.1 derivation of expressions for cooling load 

A cold storage oan be of any shape like cubical 
or cuboidal. In order to generalise the derivation it 

is assumed that the cold storage is cuboidal in shape 
as shown in fig. 2.1. 


that 


Det the volume of the cold storage be V such 


v - a . b . h 


(2-1) 


wnere a,b and H are the length (along east-west), width 

and height respectively. The sink below the ground level 
is denoted by s* 


The total heat-load is given by 

t Q„ + Q. -L Q 
1 2 w inf + rixed 


( 2 - 2 ) 


where Q, and Q 2 are heat transfer rates above and below 
the ground level respectively. is the infiltration 
heat load and Q f lxe<3 is the fixed heat load, that is 

independent of the insulation thickness, sink and the 
dimensions of the cold 


storage. 
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The heat transfer rate above the ground level 
is given bys 


1 ~ * a ' (H “ s) L2 A T + S at +S aT 

o gg 

+ U - b(H - s) CEiT E + Eat se + + s*T w ] 


OlM * 


+ U roof * a • b [2 AT R + £ aT^] 


(2-3) 


Sre Tse 7 T SS ' T sw and T sr are th e magnitudes of the 
temperature above ambient value to account for the 

radiation effects (sol-air temperature) [ 3 ] , for the 

Eastern, Southern, Western sides and the roof respec- 
tively. 


A 1 > — T »_ rp 

ss “ X SS T s 


(2~4a) 

AT - T _ rp 

ss ” 5 S ~E 


(2 -4b) 

AT ss T _ t> 

sw 1 sw T w 


(2— 4c) 

AT — T _ rn 

SR “ X SR T R 


( 2— 4d) 

The overall heat transfer coefficients 

are given by 

U ■— - rill 1 _ i 

k 

(2-5) 

rl 1 {■’ IT — f — t: — = 

% + h + | + 2 (V k l)] 

1 o t=l 

and. 

t +t~ 
e 

u rn . - J 

k 

( 2—6 ) 

roof 4 - rs . 4 _ ***^ rt *^ ***** rr 

r 1 1 S n*t 

Ck I + % + + mli'rd 

° m 



where and t' are equivalent thickness for the 
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n— layered walls and n* layered roof, respectively. 

The thickness of the main insulation and the 
conductivity are t and k„ respectively. 

Therefore, 


Q i = '(wt (H ~ s) [a cs AT n + s A Tg + s A Tgs . 

sat sn] + bl s at e + s at se tuy + 

+ TtT+tT • a . b [Eat r + E« a ] (2-7) 

jR 0 

The heat load for the underground structure is calculated 
as : 


0 2 = U D 2 (a + b) S E at i 


(2-3) 


where is given by 


U D = 


t„+t 


r 1 1 n t. t ~l , 'e 

[h7 r " + + ( k^ + ir 3 

(2-9) 


The magnitude of b." for the underground structure is 
assigned to be 2 k J m^h °k based on practical conside- 
rations [ 13 ] . 


The quantity AT^ is calculated from the under- 
ground temperature data [ 19] . 
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Thus 0 2 = £-~p£7> 2 (a + b) S Sat (2-10) 

1 e 

The infiltration heat load is given by 
x a « b (H-s) ^ch v,/, 

Q inf ~ v"' " “ 24 ' s ^ h a (2-11) 

where i 3 the number of air changes per day given by 

the empirical relation 

D c h = 110 * 03 - 54.3906 In (V) + 11.4-26 [ln(V)] 2 
- 11.32 [ln(V)] ? +0.043703 [ln(V)] 4 


( 2 - 12 ) 

^a / h i are the enthalpies of the ambient air and inside 
air respectively. 


The design heat load is then given by 


^des ~ + ^i n f + ^2 

Over the peak load hours i.e. 10 a.m. to 3 p.m. 


(2-13) 


3pm. 


. r v 3pm* 

Q des = 1,1 L Tt+T") (H “ s) * a ( S AT n +SAT 


3prc, 

+ S 

10 at 

3pm. 

: 2 A 

10am. 


+ TtL+ET a 

K. ■ <3 


10am, 

3pm. 


10am. 


» + 

E A T_ e ) 

SN 

10am. SS 

3pm. 

3pm. 

+ ^ A t se 

+' 2 A % + 

10am. 

10am. 

3pm. 

3pm. 

b . [ E a 

T + S a T 

10am. 

K 10am. 


3pm. 

S 

10am ( 


'SST 


SR- 


eqn. contd. 
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1 c . 3pm. 

+ ~HF Tp>) 2 (a + b) S Eat. 
(t l + V 10am. 1 


a . b (H-s) n ch 3 P m - 

+ ... . '"24 ^ < h a - h i>J < 2 -M) 

a 10 am. 

^fixed has not been considered here because it is inde- 
pendent of t / t R / t 1 , a, b, H and s. 

The factor 1.1 has been incorporated in the abov 
equation to account for the heat transfer through the 
floor [ 20 ]. 

The tonnage of the refrigerating system having 
SP as a safety factor is given by 


SP Q, 

n _ _ _Jes 

R ~ 127600 x~5 


(2-15) 


The actual heat load is given by 


W " ^ [ TETO' <H - S) la <i * T n fi * T S 

e 24hrs. 24hrs. 


+ E ^ T cw t E AT _) 

SN 24hrs. 


24hrs 


+ b ( 2 A T + 2 aT + 2 aT 

24hrs. * 24hrs. SE 24hrs. 


+ E aT ) i 
24hrs. 


a . b { 2 AT + EaT j 
ir R + V 24hrs. R 24hrs. SR 


+ a— Tirr>) 2 (a+b) S 2 a m 

Ct l + *0 ] 24hrs. 1 


eqn . contd . 
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+ P ch 

Va ^ 2 Jh <h ®' h i ,] 

2.1.2 HEAT TRANSFER 3Y qnr^ ^ 

riY S0LA R Radiation 

In the problem of heat transfer «, 

“* -of, the two princioal f t TOllS 

' cipal factors of the evt^ , 

thermal environment are the out do < 

, c3oor a * r temperature 

solar radiation intensitv Both 
erratic fi *. 3 < °th are subjected to 

ratxc fluctuations of w h lc h the former folio 
tiaii, 7 , ~°Ho v/s essen- 

ally a periodic variation on clear a 

Present enalvsis a . ^ In «» 

- slS/ a instant internal +->, 

and periodic variations f ’ * envlroI ™ent 

solar radi f ^ ^ and 

radiation intensity are considered. 

Figure 2.2 shows a diurn ->i 

door air t Nation in out 

door axr temperature for Kanpur citv drr- 

hot day in June. It t * tYPiCally 

t can be obse^^v^r? +-v> 4 - * 

rature usually occurs Just before _ ^ 

temperature usually 0 * lle 

7 CUrS SOne 1 to 2 hours after the 
solar noon. e 

Figure 2.3 exhibits 

mbits the variation in direct 

solar radiation over the dav. It 

,, X1Z can oe seen that the 

direct solar radiation reaches a m ■ 

and ■ 1 "aximum at solar noon 

and exists onlv Airim 

xy Gur; i-ng sunrise to 


sunset 
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2.1.3 INCIDENT SOLAR RADIATION 

The total radiation, i fc reaching a terrestrial 
surface is the sum of the direct solar radiation i D , the 
diffuse sky radiation l d , and the solar radiation ref- ' 
lected rrom the surroundincf surface I 

r* 

The intensity of the direct component is the 
product of the direct normal irradiation i^ , and the 
cosine of the angle of incidence 6, between the incoming 
solar rays and a line normal to the surface. Pig. 2.3. 

Thus, 

I t = I DN 003 6 + I d + T r W//rn " (2-17) 

In the present analysis the reflected radiation 
is neglected. 

2. 1.3.1 DIRECT NORMAL SOLAR INTENSITY 

At the earth's surface on a clear day l_. T is 

” ON 

represented by, 

I DN = -eapTB/slnT) (2-18) 

where A is the apparent solar irradiation at air mass 
equal to zero and 'B' is the atmospheric extinction 
coefficient tabulated in [7 ]. 
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2. 1.3. 2 DIFFUSE SOLAR RADIATION 

The diffuse radiation falling on any surface 
consists of radiation from the sky and part of the 
reflected solar radiation from adjacent surfaces, parti- 
cularly the ground lying south of the surface in ques- 
tion. A simplified general relation for the diffuse 
solar radiation I from a clear sky that falls on any 
surface is given approximately by the following equation 

X DS = C I DN F SS (2-19) 

where /C' the diffuse radiation factor, is tabulated 
in [ 7 ] and F sc , is the angle factor between the surface 
and the sky. F^^ is 0.5 for vertical surfaces add 1.0 
for horizontal surfaces. 

2.1.4 DETERMINATION OF SOLAR ANGLES 

The sun's position in the sky can be most con- 
veniently expressed in terms of solar altitude £ , above 
the horizontal, and the solar azimuth 0 , measured from 
the south. These angles /3 and 0 depend on the latitude 
of the place L ; the solar declination § , which is a 
function of the date and hour angle h. 

The following equations relate £ and 0 to the 
three angles mentioned above 

sin 3 = cos L cos 6 cos h + sin L sin 6 

. . . .( 2 - 20 ) 
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cos 0 = (sin 3 sin L - sin 6 )/(cos 3 cos L) 

( 2 - 21 ) 

The hour angle is zero ‘degree at solar noon and increases 
by (360/24)° every hour. 

The value of solar declination as a function of 
date is tabulated in [ 7 ]. 

The solar position angles and incident angles 
for horizontal and vertical surfaces are shown in fig. 
2.4, where line OQ leads to the sun, the north— south 
line is SOM and the east-west line is EOW, line OV is 
perpendicular to the horizontal plane in which the solar 
azimuth, angle HOS(0) , and the surface azimuth, angle 
POS('Jf) are located. The surface solar azimuth (r) is 
the angle HOP . The angle of incidence © for any surface 
is defined as the angle between the incoming solar rays 
and a line normal to that surface. For the horizontal 
surface in fig. 2.4, the incident angle 9^ is QOV; for 
the vertical surface the incident angle © v is QOP . 

For any surface incident angle is related 
to .6 , r and the tilt angle of the surface v by 

cos © = cos 3 cos r sin v + sin 3 cos v 

... ( 2 - 22 ) 

when the surface is horizontal v = 0, and 


cos © H = sin 3 


(2-23) 
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For a vertical surface, v = 90° and 

cos © v = cos 0 cos r (2-24) 

2.1.5 OUTSIDE HEAT TRANSFER COEFFICIENT 

The evaluation of summer thermal load for a 
given structure requires knowledge of each of the thermal 
resistance interposed between inside contained atmosphere 
and outside air. One such resistance is that of the film 
adjacent to the outside surface. 

By common usage this outside combined (radiation 
plus convection) film coefficient is taken as 122.7 kJ/m — 
for summer [ 6 ]. The above standard film coefficients 
are reasonable average values and their use permits tabu- 
lation of overall coefficients of heat transmission for 
a multitude of types of construction. However it should 
be recognised that, in particular localities there is a 
marked variation from the standard value* This section 
presents a simple tabular method for obtaining a design 
value of the outside heat transfer coefficient as a func- 
tion of local climatological conditions, and correct- 
ing the " standard " values so that they will be more 
precisely applicable to different locations. 

li is given by the following empirical rela- 
o 

tions. 
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For very smooth, smooth, moderately rough and 
rough surfaces h Q is respectively given by. 


h 

o 

= 14.273 + 3.559 V . , 

wind 

(2.25a) 

h 

o 

= 13.42? + 3.31 V . , 

wind 

(2.25b) 

h 

o 

^ 26.55 + 5.0663 V . _ 

wind 

(2 ,25c) 

h_ 

o 

s 23.64 + 6.364 V . , 

wind 

(2.25d) 


2 

in IJ/m -h-°c when v w j_ n( g is given by }cm/h, 

2.1.6 CALCULATION OF SOL-AIR TEMPERATURE 

In heat transmission calculations, it is con- 
venient to combine the effects of outdoor air tempera- 
ture and solar radiation intensity into a single quan- 
tity. The rate of heat transfer from the external 
thermal environment to the outside surface of a sunlit 
wall or roof may be written as 

a = h (T-T ) + a I (2-26) 

X) o o w,o s 

where T is the teraoerature of the outdoor air, T 

o - w,o 

is the temperature of the outside wall surface, f% is 

s 

the solar absorptivity of the wall, and I is the combined 
incidence of solar radiation (direct and diffuse) upon 
the surface. The rate of heat transfer may be expre- 
ssed as 

a = h (T - T ) 
t) o e w,o 


( 2-27 ) 
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Prom Eqns. (2-26) and {2-21), one gets 

T e = T o + a 1 - /h 0 (2-23) 

This fictitious temperature T g is called sol-air 
temperature. 

The variations of sol-air temperature with time 
may be assumed to be repetitive for successive 24 hour 
cycles. Any of the curve may be mathematically e>q>ressed 
in terms of a Fourier series. Thus if T & = f(e), where 
9 is the number of hours measured from midnight solar 
time, we have 

T e = T e,m + M i 003 “i e + iq 1 sin “j© + cos w 2 © 

+ N 2 sin W2 © +..... (2.29) 

where the coefficients T , M and N are given bv 

s / in n n 


, 24 

~ r 

" 24 0 

T 

e 

de 


(2-30) 

_i ? 4 

“^o 

T 

e 

cos 

w e ae 

n 

(2-31) 

, 24 

= 12 ' 
o 

T 

e 

sin 

wed© 

n 

(2-32) 


where oj- - t/12 radians per hour and ,,, = nw„ . 

1 w n 1 

2.1.7 INFILTRATION LOAD 

Infiltration load is due to the air that leaks 
into the conditioned space from the surroundings are 
viceversa. The exact estimation of the leakage of air 
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through apertures is extremely difficult. However,, 
approximate values have been tabulated by ASHRAE [ 7 ]. 
Similarly the leakage of air due to door openings and 
through shutters is also based on practical experience 
in terms of the room sizes. The volume of infiltrated 
air is related with the volume of the room for different 
usages - average usage, long usage and heavy usage. The 
infiltration heat load is already described by equation 
( 2 - 11 ). 


2.1.3 HEAT TRANSFER THROUGH THE STRUCTURE 

The heat transferred through a multilayered 
structure of a system having N layers from the outside 
environment to the inside is given by 


Q = UA (T - T. } (2—33) 

0 . 1 

where U is the overall heat transfer coefficient given 
by 


i N t. 

Jr. _ i J. y i. . JL. , 4_ b 

u - h + k7 + h + k 

1 1-1 1 o 


(2-34) 


where t is the insulation thickness and k is the con- 
ductivity of the insulation material. T and T. are 

a l 

the ambient and inside temperatures, respectively. A 
is the area of the wall through which heat is being 


transferred 



23 


/ the inside heat transfer coefficient is a 
function of the temperature drop across the film formed 
ov^r the wall surface, is given by the following 

empirical relations [6] 

h. = 5.593 (AT) 0 * 25 ( 2— 35a) 

for floor h i = 6.0276 (aT) 1//3 (2- 35b) 

for the ceiling 

where tu is in kJ/m 2 ~h-°c and 

AT = <T wi - T.) (2-36) 

The above equation for hu cannot be solved directly/ 
since in .such case the film coefficient is a function of 
at across the film/ and at is itself a function of the 
film coefficient. An iterative trial and error procedure 
is therefore necessary. 

2.2 ECONOMIC MODELS 

The following four economic models are generally 
considered in the analysis of any investment policy. Any 
of them can be used to indicate the most favourable 
choice amongst several alternatives. Some of the merits 
and demerits of the same are briefly discussed. 
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2.2.1 PRESENT WORTH METHOD 

In this method all the costs are transferred 
into the present worth. In order to determine the pre- 
sent worth of a future sum , the interest rate is 
established first. Then present worth is the value of 
a sum of money at the present time that, with compounded 
interest, will have a specified value at a certain time 
in future. Thus 

. P _ = (PWF) S (2-37) 

( 1+r) n 

where P = principal 

r = rate of interest per period 
S = total amount to be repaid at future time 
n - number of periods 

PWF = single paavment present worth factor. 

The series present worth is given by 

P = R [ ] = (R) (SPWF) (2-33) 

r( 1+r ) n 

where P = original amount 

R - amount withdrawn at the end of first and 
subsequent periods and 
SPWF s= series present worth factor. 

Although the present worth method is most widely used, 
it presents some difficulty if the life times of two 
possible investments are different. 
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2.2.2 ANNUAL COST METHOD 

It transforms all non-annual costs to an annual 
basis. . In this method also the rate of Interest must be 
established first. 

R = P [ ] = P (CRP) (2-39) 

(1 + r) n ~l 

where CRP is the annual recovering factor. 

The annual cost method has two advantages over the pre- 
sent worth method. One is that there is no complication 
introduced when the prospective investments have diffe- 
rent livss. The second is that it is more natural for 
most people to think in terms of an annual cost than in 
terms of the present worth. 

2.2.3 RATE OP RETURN METHOD 

In this method no assumption is made regarding 
the interest rate. Instead the money for the investment 
is considered to be in hand/ and the rate of return is 
calculated as though it were an interest rate received 
by making an outside investment. Here we find r such that 

P (CRP) - salvage value (SPP) = Net income 

where SFP = Sinking fund factor = [ r — -] (2-40) 

(l+r) n -l 



31 


2.2.4 BREAK EVEN POINT 

This method assumes that money is borrowed at a 
specific rate of interest and that the loan is paid off 
as rapidly as possible with no profits extracted. The 
break even point is defined as the time where the loan 
is paid off and profits begin following to the investor. 
In contrast to the rate of return method where the inte- 
rest rate was unknown, the life is unknown in this method. 
Thus one finds n such that 

P(CRP) — salvage value (SFE) = Net income 

(2-41) 

2.2.5 ANNUAL COST X-fETH SIMPLE INTEREST 

A simple economic model of annual cost with 
simple interest instead of compound interest can also be 
considered. Here 

r _ (2.42) 

In the present work, the present worth method is consi- 
dered to decide about the investment policy for the 
optimization problem. 

2.2.6 COST ANALYSIS 

The various individual costs which contribute 
to tlie total cost in terms of present worth follows: 
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(l) CD ST OP INSULATION FOR WHOLE STRUCTURE: 

c^[2(H-s) (a+b) t' + 2s(a+b) t^ + a b (t^ + t^) ] 


C 1 = 


... (2-43) 

where l s the cost of insulation per unit volume and 
L is the life of machinery and insulation. 


(2) COST OF POWER: 


C 2 “ 


C Q . -> x 365 0 r- 

__eff actual ^ _____ 3.5 

12/600 ~ * COP 


where 


'eff 


1 

L 


e i 


L 

S ■ 
i=l (1 + R) 


i-1 


(2-44) 


(2-45) 


. th 


c is the cost of electricity for the i vear 
e i 

(Appendix - B). c e ff is the effective cost of electri- 
city (Rs/KWH). F is a factor for the actual power con- 
sumption of the system per year. R is the rate of inte- 
rest. 


(3) COST OF REFRIGERATING MACHINERY: 

c _ C 3 * 561 ^ desig n - (2-46) 

3 5 X 12/600 x L 

where c_ is the cost of refrigerating machinery per ton 
capacity. SF is the safety factor for refrigerating 
machinery. 
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(4) MAINTENAN CE COST: 


C 


4 


( 0 . 10 ) 


„ r (i + r) l - l i 

c |_ — -• j 

D R(1 + R) L_1 


(2-47) 


maintenance cost has been considered as 10£ of the equip- 
ment cost [ 20 ]. 


(5) DIGGING AND ADDITIONAL UNDERGROUND CONSTRUCTION 
COST: 


Cj- = c- (s . a . b) (2-48) 

5 D 

where c- is the digging and additional cost for the 
structure per unit volume. 

The total cost per year is then given by 

-EC. (2-49) 

i=l 1 




GIAPTER_ _3 
OPTIMIZATION 

3.1 TERMINOLOGY AND STATEMENT OF THE PROBLEM 

Optimisation is the process of determining the 
best results under the given circumstances. 

The optimisation problem is stated as follows s 

Find X ={ x 1 , x 2 , x n j (3-1) 

which minimises f(X) subject to the constraints 

gj(X) < 0 , j = 1,2/ ... ,m { 3— 2a) 

and 1 .(X) “ 0 ' J = m+2/.../P (3~2b) 

where X is an n-dimensional vector called the design 
vector., f (X) is called the objective function and g^(X) 
and lj(X) are the inequality and the equality constraints/ 
respectively. 

The design vector describes a set of quantities 
of an engineering system which are viewed as variables 
during the design process. The design variables are 
collectively represented as a design vector 

• X =£ x 1 * x 2 x n i 


( 3-3) 
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In practical problems, the design variables 
cannot be chosen arbitrarily, instead they have to satisfy 
certain conditions and requirements. The restrictions 
that are placed on the design variables in order to 
achieve an acceptable design are called design constraints* 

It can be seen that in most engineering design 
procedures there exists more than one acceptable design. 
Since the purpose of optimization is to choose the best 
out of the many designs, a criterion has to be chosen for 
comparing the different alternate acceptable designs and 
selecting the best one. This criterion with respect to 
which the design is optimized, when expressed as a func- 
tion of the design variables is known as Objective 
function' . 

In the present study the objective function is 
the total average annual cost of the cold storage which 
comprises the initial and running costs based on the 
present worth method. The design variables are the 
thickness of insulation on the walls above the ground 
level and on the ceiling, the sink (depth below ground 
level of the building) and the length and breadth of the 
cold storage structure. The insulation thickness for 
that part of the walls which is below ground level is 
taken as some fraction of the optimum insulation thickness 
for that part above the ground. The height of the cold 
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storage structure is a function of the length and breadth 
for a given volume. 

Thus our object is to determine the insulation 
thickness/ sink and the dimensions of the structure in 
such a way that the average annual total cost is minimum. 

The constraints on insulation thickness are 
characterised by the form in which the thermal insulation 
panels are available in the market. Some insulants like 
cdrkboard and polystyrene foam are available in sheets 
of standard thickness while other insulant like glass- 
wool is available in loose form and has to be compacted 
into sheets or slabs. 

The constraint on the sink is governed by the 
fact that the sink has to be positive and that it cannot 
exceed the height of the structure. The extreme condi- 
tions are the whole building being above or below ground— 
1 evel . 

The constraints on the dimensions of the cold 
storage depend on the volume of the cold storage and on 
the fact that volume is in itself governed by the amount 
and specific volume of the stored commodity. The maximum 
length of the cold storage is limited from the point of 
view of handling and good air distribution. The height 
of the cold storage is decided from the height of each 
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pallet and the number of pallets to be stacked-* Addi- 
tional height for safety, manoeuvring and air distribu- 
tion is also accounted. 

With these considerations the height of structure 
may vary from 6 to 15 m and the width varies from 15 to 
25 m. 


3.2 OBJECTIVE FUNCTION 

In the present stxady the objective function is 
the average annuel cost for the entire life span of the 
insulation and machinery in terms of the present value. 

The average annual cost comprises both the initial and 
running co sts . 

The optimization procedure is independent of the 
fixed heat load and hence the fixed costs. Fixed load 
takes into account the heat load from the stored commodity, 
heat generated inside the conditioned space and venti- 
1' ation heat load. 

It can be seen that the optimization variables 
like insulation thickness, sink or dimensions of the room 
do not influence the fixed load and hence the fixed cost. 

In this regard the objective function does not comprise 
the fixed cost. 
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The various individual costs incurred (in terms 
of the present value) over the entire life span of L years 
have been divided by L to compute the objective function. 

The initial cost comprises the cost of insula- 
tion, the cost of machinery for refrigeration, the cost 
of digging and underground constructional cost. The 
running cost comprises the cost of power (electricity) 
and the cost of maintenance. 


The cost of refrigerating machinery is a func- 
tion of the tonnage for the required system. The tonnage 
of the system is given by 
SP . 0.. . 

TR ~ tons (3-4) 

12,600 x 5 

where Q, is the design heat load for five hours from 

design 

lO AM to 3 PM. 


The cost of po^'^■er is a function of the actual 
heat load (^ ac ^. ua ^) for 24 hours and the coefficient of 
performance of the system. 


The maintenance cost is assumed to be 10/ of the 
cost of refrigerating machinery [ 20 ]. 

The cost of insulation and digging are functions 
of the dimensions of the cold storage. 

The objective function is the average annual 
cost being expressed as the present worth of all the 
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expenditures . 


f (X) = C. .. . , + C 

initial running 

= (C 1 + -+ C 5 ) + (c 9 + C 4 ) (3-a) 

The functional dependence of the objective function is 
expressed as 

Average Annual Cost = f (t,t_ , twaAH, s) 

K 1 

+ g( constants) (3-6) 


Our object is to minimize this average annual cost by 
varying t A t^,a,b and s subject to constraints. 


3.3 OPTIMIZATION TECHNIQUE 

An optimization problem in which there is no 
restriction on the choice of the design vector, is called 
an unconstrained optimization problem. Any restriction 
on the design vector constitutes what is called the 
constrained optimization problem. Most of the problems 
in engineering applications are constrained, but some 
of the most powerful and convenient methods of solving 
constrained problems involve the conversion of the cons- 
trained problem to one of unconstrained minimization and 
this technique has been used in the present work. 

3.3.1 THE UNCONSTRAINED MINIMIZATION TECHNIQUE 

The multidimensional unconstrained minimization 
technique being used in the present work is the Powell's 
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conjugate direction method. 

Detailed description of this method is available 
in any standard text on optimization [21 ]. 

In this method at each point in the function 
space under consideration/ there is a preferred direc- 
tion along which the values of the design variables are 
changed systematically by a well defined scheme/ which 
leads to the minimum. Thus starting from the initial 
guess one arrives at the minimum through the sequence 
of improved approximations/ each derived from the pre- 
vious approximation. The Powell's method is a non- 
gradient method and it has the property of quadratic 
convergence. Without going into any mathematical detail 
a brief description of the algorithm is presented here. 

3.3.2 POWELL'S METHOD 

It is an extension of the idea of pattern move. 
For an intuitive understanding, Fox [ 22 ] describes the 
method as follows t Given that the function has been 
minimised once in each of the co-ordinate directions and 
then in the associated pattern direction/ discard one of 
the co-ordinate directions in favour of the pattern 
direction for inclusion in the next m minimizations, 
since this is likely to be a better direction than the 
discarded one. After the next cycle of minimizations. 



41 


generate a new pattern direction and again replace one 
of the co-ordinate directions. Pig. 3.1 explains the 
algorithm. The basic method has a tendency to- choose 
nearly dependant direction in ill-conditioned problems 
and the method may fail to converge to the actual minimum. 
One simple remedy is to reset the directions to the origi- 
nal co-ordinate vectors periodically and whenever there 
is some indication that the directions are no longer 
productive Powell has recommended a more effective proce- 
dure to overcome this difficulty [ 23 ]. The search is 
terminated when the relative change in the function value 
and in all the design variables/ between two consecutive 
cycles of minimization/ is less than the desired accuracy. 


3.3.3 QUADRATIC INTERPOLATION 


In this technique the function F(oc) is approxi- 
mated by a function H(a) which has an easily determinable 
minimum point. H(a) is expressed as 

H(a) = a + boc +• c« 2 (3-7) 

the minimum of which occurs at 


dH 

dec 


b + 2coc = 0 


(3-3) 


0* = - b/2c 


(3-9) 


The constants a, b and c for the approximating quadratic 
can be determined by sampling the function at three 
different a values, cc^ , and and solving -the 
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equations 

Fj = a + ba 1 + ca^ 2 (3-10) 

f 2 = a + ha 2 + ca 2 ( 3 - 11 ) 

and = a + bcr.^ + coo ? 2 (3-12) 

where F. , F 0 and F stand for F(cc.,), F(a 0 ) and F(a_) 
r e sp ec ti v ely . 

Solving the above three equations, the values 
of a, b, c and <X* are obtained as 


T( W Vj) - 


VJ ~ 

[(oc t -a 2 ) (a 2 -a^) (a^-a^) ] 



* • # 

(3-13) 

b = 

(F 2 ~F 1 )/ (ff^-^) - c(a 2 +a ) 

(3-14) 

o ~ 

f rV“i 2 

(3-15) 

a 

~b 

2c 

(3-16) 


Choosing a, r and 2r for o:^ , oc^ and respectively, 
where r is the preselected trial step-length. One 
obtains 



a 

-TP 

"i 

(3-17) 


b 

= (4F 2 - 3F 1 - F,) / 2r 

(3-13) 


c 

= (F„ + F t - 2F 2 )/ 2r 2 

(3-19) 

and 

a 

(4F 2 - 3F t - F 3 ) 

- xW 2 ' j -'Tf^“2f“) 

(3-20) 
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For a to correspond a minimum for H ( a ) , 
it must satisfy 



> 0 

ascc* 


(3-21) 


Since H is a quadratic, this requires c > 0 
or (F^ + F 1 ) > 2F 0 . 

The point oc* is considered to be a sufficiently 
good approximation to the minimum of F(a) if 

H(a*) - F(a*) 

~ — < s (3-22) 

F(a*) 


where e is a small number to be specified depending on 
the accuracy desired. Fig. 3.2 explains the algorithm. 


3.3.4 THE PENALTY FUNCTION METHOD 

The basic function of the interior penalty 
function method is to convert the original constrained 
problem into one of unconstrained minimization by blend- 
ing the constraints into a composite function and making 
it possible to ignore them at the minimization stage* 

In this method the numerical solutions are sought by 
solving a sequence of unconstrained minimization problems. 

The advantage of penalty function method is 
that ^powerful , well studied and reliable algorithms for, 
the unconstrained minimization of arbitrary functions 
can be used in this method. 
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The penalty function formulations for inequality 
constrained problems can be divided into two categories; 
interior and exterior. In the interior formulation, the 
unconstrained minima lie in the feasible region and 
converge to the solution as a special parameter is varied. 
In the exterior formulations, they lie in the infeasible 
region and converge to the solution from the outside. 

The advantage of the interior penalty function 
method is that, given an initial acceptable, it produces 
an improving sequence of acceptable decisions. Moreover, 
we approach the constraints in such a way that they become 
critical only near the end. This is a desirable feature 
in engineering design because one is free to choose a 
suboptimal but less critical design if required. 

After studying' the relative advantages and dis- 
advantages of the available methods, the interior penalty 
function is used in the present work. 

In the interior penalty function method a new 
function (0 function) is constructed by augumenting a 
penalty term to the objective function. The 0 function 
is defined as 

m 

0 (D, rp - p(d) - r* J E 1 jTfcy 

• 3 

where F is to be minimized over all D satisfying 


(3-23) 
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g j (D) < 0 , j = (3-24) 

The penalty parameter r is made successively smaller in 
order to obtain the constrained minimum of F. 

The flow diagram of fig. 3.3 explains the 
algorithm. t 

The interior penalty function method requires 
a feasible starting point for the search towards the 
optimum point. In cases where finding a feasible start- 
ing point is difficult, the penalty function method 
itself can be used for finding it, but the procedure 
is time-consuming. 



CHAPTER - 4 

RESULTS AND DISCUSS ICN _S 

4.1 COMPUTATIONAL METHOD 

The objective function formulated in Chapter 2 
has been optimized, as outlined in Chapter 3/ for various 
inside temperatures and insulation material. 

An Ammonia Vapour-Compression refrigeration 
system has been considered for the necessary cooling of 
the cold storage space. Life of the insulation and the 
refrigerating system have been taken as 20 years and 
an operating factor of 0.75 for power consumption has 
been considered. A condenser temperature of 42 °c and 
evaporator temperature of - 3 v c have been assumed. 

Further it has been assumed that the refrigerant is being 
superheated by 5“c and undercooled by 5°c. Under these 
conditions, the coefficient of performance of the system 
turns out to be 3.5511. Empirical relations for refri- 
gerant properties are given in Appendix-A. 

A method to compute the present worth effective 
electricity charges by forecasting future electricity 
costs is discussed in Appendix— B [ 2 -1 j . 
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A computer programme was made, based on mathe- 
matical formulation, having a generalized approach and 
realistic methods. The formulation deals with two cases 

(i) Determination of optimum insulation requirement 
for a given size of the cold storage and 

(ii) Determination of optimum insulation thichness 
and cold storage dimensions for a given volume 
of the cold storage soace. 

Though the programme is quite general, typical 
results are displayed for the following data? 

1) Volume of the cold storage = 6000m'". 

2) Storage conditions = 0°c and 35/ relative humidity. 

3) The upper bounds on the length, breadth and height 
being 30 , 30 and 10 m respectively. 

4) Two t^’pes of insulating material namely polystyrene 
foam and thermocole. 

The properties and costs of these insulants are 


Insulation 


Thermal conductivity Cost per un^t 
kJ/(m h v c) volume Rs/rtT 


Thermo cole 0.175 39 1120.00 

Polystyrene foam 0.150 2000.00 


5 


\ 


Cost of refrigerating machinery for Ammonia is 
Rs . 5000/ton . 
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6) Environmental temperatures and solar radiation data 
[Appendix-c] . 

4.2 GRAPHICAL PRESENTATION OF RESULTS 

Figure 4.1 shows the variation in ambient tempe- 
rature and the sol-air temperatures over the different 
walls and the roof. 

Results arc presented for conditions prevailing 
on a typical day in June from 6 AM to 6 PM. n Exact time 
of sunrise and sunset are not made use of because, -solar 
radiation at sunset and sunrise was not available". From 
this plot it is seen that sol-air temperature for a 
horizontal surface (roof) is the highest and is maximum 
at noon, corresponding to the maximum incident normal 
solar radiation. From this plot it can be inferred that 
heat transfer through the west facing and east facing 
walls is more prominent than that for north and south 
facing walls . Thus the area of eastern and western walls 
should be less than that of northern and southern walls. 

The variation in design heat load with sink for 
a typical cold storage having an inside temperature of 
0°c, exhibits significant reduction in heat load with 
sink. Fig. 4.2. The infilt ration heat load also decreases 
with sihk / Fig.4.3. Figure 4.4 shows the nor cent age reduc- 
tion in design cooling load with sink. Percentage 
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reduction in cooling load is given by 

/ reduction ~ 100 [ Q. . ^ , . . 

des/ without sink 

£ ^des, without sink] 


Q i 

des ,with sinkj/ 

( 4 . 1 ) 


The insulation thickness pertains to that part of the 
walls which are above the ground level. Insulation 


thickness for the ceiling is assumed to be 1.5 times that 


above ground level and thickness for walls below ground 
level as 0.75 times that above ground level. This plot 
gives an idea of the saving in cooling load with sink. 


Figure 4.5 shows the specimen variation in indi- 
vidual costs with sink. Among all the individual costs, 
the cost of power is highest. Thus the total cost can 
be controlled effectively by reduced power cost. This 
can be achieved by sinking the structure below ground as 
the cost of power reduces significantly with increase in 
sink. It is interesting to note that the cost of power 
at a sink of 9 m is only about 53# of that without any 
sink for thermocole.- This is due to the fact that the 
cost of power and the cost of refrigerating machinery 
are dependant on the actual heat load which decreases 
with increase in sink. The cost of maintenance closely 
follows the cost of refrigerating machinery. 

The only cost that increases with increase in 
sink is the digging and underground constructional cost. 
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Low cost, of digging and underground construction per unit 
volume, therefore cannot overcome the effects of the 
other costs to determine optimum sink. At low unit cost 
of digging and underground construction minimum total 
cost occurs at maximum sink. High digging costs, however, 
result in intermediate values of sink lying between the 
upper and lower bounds. 

Figures 4.6 and 4.7 shox-7 the variation in total 
cost/year with insulation thickness for polystyrene foam 
insulation and thermocole. These plots are for chosen 
values of the dimensions which are not optimum values. 
Variations in different sinks are also shows in these 
figures. The optimum thickness of insulation can be 
easily read from these figures corresponding to the 
minimum total cost/year. 

It can he observed from these plots that optimum 
insulation thickness decreases with increase in sink. 
Further, the optimum insulation thickness values for 
polystyrene insulation is considerably less than that for 
thermocole for the respective sinks. This can be attri- 
buted to the fact that the thermal conductivity of the 
polystyrene foam insulation is less than that for ther- 
mocole. The total cost corresponding to the optimum 
insulation thickness for a particular sink is found to 
be lesser for thermo cole insulation than polystyrene. 



51 


Figure 4.9 shows the variation in total cost 
with insulation thickness for polystyrene insulation. 

Here the total cost is based on the suggested design 
approach (consideration of hourly data for temperature 
radiation). Further the plot is drawn for assumed optimum 
dimension, i.e. length and breadth. For a given volume 
and fixed height, it was seen that total cost is minimum 
for equal values of length and breadth. This can be 
explained by the fact that for a given volume the cubical 
shape renders minimum surface area. Minimized surface 
area results in lesser cost of insulation and lesser 
heat transfer through ceiling. 

Figure 4.9 shows the variation in total cost 
with insulation thickness for thermocole insulation, 
based on suggested design approach for optimum dimensions. 
The saving in total cost for the optimum insulation thick- 
ness, for a given sink can be determined by comparing this 
clot with Fig. 4.7, where the dimensions are not optimized. 
As an illustration a saving of approximately 1/ in total 
cost was obtained for optimum insulation thickness values 
for tho case of zero sink (when the whole structure is 
above erround ) . This sav?.ng was seen to decrease with 
increase in sink. 

A saving in total cost can bo obtained if 
thermocole is used instead of polystyrene 4‘S^lation for 

M* M. rt 
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optimized dimensions as can be seen by comparing Pig. 4.3 
and Fig. 4.9. As an example a saving of 13 / was obtained 
for the optimum values corresponding to an intermediate 
sink of 5m. 

Figure 4.10 depicts the variation in total cost 
with insulation thickness for polystyrene foam based on 
standard design conditions (peak value design) . It can 
bo ascertained that the suggested design approach is 
more realistic, by comparing this plot with that of 
Fig. 4.3 which is based on the suggested design approach . 

Other parameter remaining the same, a saving of 23/ is 
obtained in cost corresponding to the optimum insulation 
thickness at a sink of 0 m (whole structure being above 
ground level). Further the saving in cost increases 
with increase in sink. It is seen that for a sink of 
9 m the saving is 42/. 

Similar savings in total cost for thermocole . 
insulation can be ascertained from Fig. 4.11 in compari- 
son with Fig. 4.9. Here a saving of 33/ is obtained 
corresponding to the optimum insulation thickness at a 
sink of 0 m, based on suggested design approach. A sav- 
ing of 46/ is obtained for a higher sink of 9 m. 


Thus it can inferred that the standard design 
procedure is more severe in the case of thormocole as 
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compared to polystyrene, as far as saving in total cost 
is concerned. The variation in total cost with sink for 
foam insulation when the size of the cold 
storage is fixed can be studied, from Fig. 4.12. 

It is seen that at low insulation thicknesses 
the optimum sink lies on or very near the' upper bound , 
for sink^ which is 9 m. However, the optimum point moves 
away from the upper bound at higher values of insulation 
thickness. This can be attributed to the fact that insu- 
lation cost overcomes the effect of digging cost on the 
total cost. 

. The variation in total cost with sink for 

thermocole insulation, for the case of the cold storage 
size being fixed is shown in Fig. 4.13. It can be seen 
that optimum sink values, at higher insulation thickness, 
are less for thermocole insulation than that for poly- 
styrene insulation (Figs. 4.13 and 4.12). 

The variations in total cost with sink for 
polystyrene and thermocole insulations, based on the 
suggested design approach for assumed optimum length and 
breadth is shown in Figs. 4.14 and 4.15, respectively . 

The optimum length and breadth has been assumed on the 
basis that for a given, volume the cubical shape renders 
minimum surface area. 
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The variations in total cost with sink for 
polystyrene and thermocole insulations on the basis of 
standard design conditions is shown in Pigs. 4*16 and 
4.17, respectively.. The total cost corresponding to 
any insulation thickness for the optimum sink is always 
higher for polystyrene than thermocole. By comparing 
Pigs. 4.14 and 4.16 it can be seen that a cost saving of 
34^f can be achieved by the suggested design approach, 
in the case of polystyrene insulation for an insulation 
thickness of 8 cm. Similarly, a cost saving of 38 / can 
be achieved for thermocole insulation for the same insu- 
lation thickness. 

4.3 TABULATED RESULTS 

Table 4.1 presents the variation in optimum 
insulation with inside temperature and sink for thermo- 
cole insulation. It is observed that for the same inside 
temperature the design heat load decreases with increase 
in sink. But this is not true in the case of total cost 
variation with sink for all inside temperatures. At some 
inside temperatures (for cooling) where the cost of power 
and refrigerating machinery is relatively low, the cost 

of digging and additional underground constructional 
may overcome the cost of power 
Cost/in increasing the total cost at some sink. In such 

case the optimum sink lies between the upper and lower 

bounds of the sink. 
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This "table also gives an idea of the variation 
in optimum insulation with sink. 

Table 4.2 shows the variation in optimum insula- 
tion with inside temperature and sink for polystyrene 
foam insulation. 

Results have been tabulated for inside tempera- 
tures of 10°c and 4°c (for cooling) and 0°c and — 4°c 
(for freezing) . It is seen that in the case of inside 
temperature being 10°c the value of the optimum insula- 
tion thickness increases after a particular sink indicat- 
ing an optimum sink between the lower and upper bound of 
0 and 9 m., respectively. However in all the other cases 
the optimum sink lies very near the upper bound. This can 
be explained by the fact that the cost of power and refri- 
gerating machinery is relatively low for higher inside 
temperature than for lower inside tempera tures . For all 
the inside temperatures considered the total cost is 
higher for polystyrene insulation than thermocole. This 
can be seen by comparing Tables 4.1 and 4.2. 

A comparative study of different orientations, 
insulants and design approaches is presented in Table 
4.3. It is seen that orientation has a small effect on 
design heat load, total cost and the optimum values. 

The optimum insulation thickness is lesser for the 
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«-S. E-M orientation, for all the cases, than NE-sw, 
NW-SE orientation. This observation helps the designer 
choose the appropriate orientation. 

The suggested design approach is more realistic 
than the currently practised standard design procedure. 

In case of thermooole a saving of 35-36/ in total cost 
acnieved for either orientation if the suggested design 
approach is used. The saving of 30-31/ is achieved in 
case of polystyrene insulation. 


It is observed that the optimum sink lies near 
the lower bound, for both the insurants and orientations 
considered, when the design approach is based on peak 
load data, because the cost of digging and additional 
underground construction is quite low in comparison to 
the other costs as they do not decrease much with increase 
in sink so as to affect the total cost. 

In case of suggested design approach, for either 
orientation and insulant, the optimum values for length 
and breadth are almost equal resulting in minimum surface 
area* But in case of standard design approach the opti- 
mum length and breadth differ because of the greater heat 
load transmitted through the walls. 

Table 4.4 indicates significant reduction in 
total cost and design heat load for the different insulants 
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and orientations if the suggested design procedure is 
used instead of the standard approach. 

Presently energy saving is being given wider 
attention and as such it may be necessary to use higher 
insulation thickness than the optimum value. In order 
to study this feasibility the insulation thickness is 
increased by 20 to 30/, Table 4*5. This renders a saving 
of 10 to 20/ saving in energy , involving only 1-7/ increase 
in total cost, for thermocole insulation. Polystyrene 
shows higher saving in energy with higher increase in 
total cost, than thermo cole, for suboptimal insulation 
thickness. 

Out of the two suboptimal cases considered 20/ 
higher insulation seems better. If cost is a desired 
criterion, thermo cole is recommended for suboptimal 
selection. However, polystyrene is preferable when 
energy conservation is the desired requirement. 



CHAPTER - 5 


CONCLUSIONS AND SUGGESTIONS 

5.1 CONCLUSIONS 

Prom the present study the following conclusions 

are arrived at: 

1. The generalized computer programme for determina- 
tion of optimum parameters based on hourly load 
calculation being more realistic as compared to 
the peak load design has been developed. 

2 . The heat load and inf ilt ration load decrease 

drastically with increase in sink, implying thereby 
a reduction in functional energy. The cost of 
power at a sink of 9 m. is only about 58/ of that 
at zero sink. Cost of refrigerating machinery is 
also reduced by a similar amount with increase in 
sink. 

3. Optimum value of insulation thickness decreases 
with increase in sink. 

4. The optimum insulation thickness of polystyrene 
insulation is less than that of thermocole for the 
same sink. But the total cost corresponding to the 
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5. 


6 . 


optimum insulation thickness for any sink is 
lesser for thermocole than for polystyrene. 


The total cost and optimum insulation thickness 
m the case of optimum dimensions is lesser than 
that With fixed sire, for any insulation. 


A saving in total cost of 30-36/ is achieved by 
the method outlined in the present work as com- 
pared to peak value design approach for either 
insulant. 


Change xn orientation of the structure does not 
result in large variations in total cost or 
optimum values of design parameters. 

Thermocole insulation is preferred over poly- 
styrene because it results in lesser total cost 
for the same functional energy. However, if the 
insulation thickness is to be limited due to 
practical considerations, polystyrene should be 
preferred as it results in lesser optimum insu- 
lation thickness. 


In case of the suggested design approach the 
optimum length and breadth are almost equal. 
However, they differ by 10-15/ in the case of 
standard design method because of the higher 
rate of heat transfer through the structure. 


9 . 
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10. A study of suboptimal design parameters reveals 
that a saving of 10 to 20/ in functional energy 
can be achieved with 1-6/ increase in total cost 
when the optimum insulation thickness is increased 
by 20-30/ for thermocole. A similar saving of 
15 to 25/ in functional energy and increase of 
2-15/ in total cost is found for polystyrene insu- 
lation with 20 to 30/ suboptimal insulation thick- 
ness. This enables the designer to choose sub— 
■optimal insulation thicknesses for the different 
insulations depending on as to which of the energy 
and cost requirements is more important for parti- 
cular applications. 

5.2 SCOPE FOR FUTURE WORK 

1) A study involving insulants other than those consi- 
dered in the present work. 

2 ) Application of other optimization techniques to 
minimize the objective function. 

3) Computation of optimum dimensions on the basis of 
the assumption that the heat transferred through 
all the walls per unit area is equal. 

4) Determination of optimum insulation thickness for 
each of the four walls. 



61 


rep_erences 


1 . 


4 . 


6 . 


6 . 


7 . 


9 . 


10 . 


11 . 


Chiub, P.c,, " Cold Storage Indus try-growth. 
Performance and Promise" All India Seminar 
on Cold Storages , 16—17 May, 1932, PP. 1 — 3 . 


Manohar Prasad, "Cold Storage Prospect and 
Analysis", Dept, of Mech. Engineering, I.I.T., 
Kanpur . 


Threlkeld, J.L., "Thermal Environmental Engineer- 
mg , Prentice Hall, Englewood Cliff, New Jersev, 
1970, PP. 312-351. 


Mackey, C.O., and Wright, L.T., Jr., "Periodic 
Heat Flow - Composite Walls or Roofs", ASHVE 
Trans., Vol . 52, 1946, PP . 28 3-296. 

Mitals, G.P., "Transfer Function Method of 
Calculating Loads, Heat Extraction and Soace 
Temperature", ASHRAE Journal, Vol. 14, Ho. 10, 

Dec. 1972, PP. 54-56. 

Kadambi, V., and Hutchinson, F.W., " Refrigeration 
Air Conditioning and Environmental Control in 
India", Prentice Hall of India Private Limited, 
1963. 

** ASHRAE Handbook of Fundamentals", 1931, Published 
by The American Society of Heating, Refrigerating 
and Air Conditioning Engineers, Inc,, New York. 

McClure, C.J.R., " Optimizing Building Energy Use", 
ASHRAE Journal, Vol. 13, No. 9, September 1971, 

PP. 64-72. 

Bonar, H.B., "Effects of Increased Energy Costs 
on Refrigeration Warehouse Construction', ASHRAE 
Journal , Vol. 16, No. 11, November 1974, PP. 53-56. 

Spielvogel, L.G.,"More Insulation can Increase 
Energy Consumption", ASHRAE Journal, Vol. 16, 

No. 1, January 1974, PP . 61-63. 

Gupta, B.D., " Automated -Optimum Design of 
Refrigerated Warehouses and Air Conditioned 
Buildings", Ph. D. dissertation. Dept, of Mech. 
Engq., I.I.T., Kanpur, August 1977. 



62 



Refrigerated Warehouse in High Rack 
S? n st fS c ^ on# 9° m Parison Cost of Insulation and 
Ex., loitation with Those of Single Storey Cold 

fS r ? g \ S ' ProgresS in Refrigeration on Science 
417-425 ° gY/ 01 * 111 ' ' IIR (Paris >' 197 3, PP. 


13 . 


14. 


15 . 


16. 


17. 


Claesson, A. and Holmqvist, S., "storage of Cold 
ana Frozen Food in Underground Mined Rock Caverns 
~. An Energy Saving Method of Storage" Refrigera- 
tion Science and Technology, 1973, Vol . 2, PP . 

31> *** .323. 


Manohar Prasad et» al., " Optimum Insulation 
Thickness and Depth of Sink for an Underground 
Cold Storage", Productivity in Cold Storage 
Industry , Organized by NPC (Chandigarh), Dec. 16-17, 

•1 AO A ~ * 


Lckmannekim, M., and Henning er, R.H., 'Computerized 
Energy Requirement Analysis” and Heating Cooling 
Load Calculation of Building", ASHRAE Journal, 

Vol. 14, No. 4, April 1972, PP . 25-33. 

Lokamanhekim, M», "Procedure for determining 
Heating and Cooling Loads for Computerized Energy 
Calculations - Algorithms for Building Heat 
Transfer Sub-routines", ASHRAE, 1971. 

William Evers, E., "Computerized Energy Analysis", 
ASHRAE Journal, Vol. 13, No. 9, Sent. 1971, PP . 
64-72. 


13. Hutchinson, F ,W. Heating Humidifying and Load 

Analysis", Ronald Press. Co* , H.Y., 1962 ,PP. 232-33. 

19. Saluja, S.N. ct. al., "investigation of Soil 
Temperature Variation", Journal of Institution 
of Engineers India, Vol. 49, No. 5. 

20. Borbely, P., " Borbelv, P., "Economic Insulation 
Thickness", XV International Congress of Refrige- 
ration, Venesia, 23-29 Sept.,.- 1979, D 1-32. 

21. Rao, S.S., "Optimization Theory and Applications", 
Wiley Eastern Limited, 1979. 

Fox, R.L., " Optimization Methods for Engineering 
Design", Addison Wesley, Reading, Mass, 1971. 


22 



63 


23. 


the e Miniii™ D of I' f E S icient M etnod for Finding 
without Cal mil +-« un ction of Several Variables 

7(4), 1964, pp ! ’ 303_ 307 ^ Vatives Con1puter Journal, 

Vaoonr T; ^S‘ Sarj ' ■ * 0ptlmura Jesign of Multistage 

E f r i9eration Astern", M.Tech • 

DeoSb^r Sli.° f MeCh * ^S- Kanpur 


24. 



64 


APPENDIX - A 

PROPERTIES OF REFRIGERANTS 
A » 1 AMMONIA R-717: 


Cp(T) = 2.643 + 0.004643[ (T-10) + 0.05015 (T-10) 2 ]for T< 10 


Cp(T) = 2 . 14375+0 .00 3969 16 [ (T+45 )+0 .00006 (T+45 ) 2 ] for T>10 

h f (T) = 180.835 + 462.25 (T/100) + 23.7163 (T/100) 2 + 

13.8 36 (T/100) 3 + 5.36214 (T/100) 4 

h (T) = 1443.36 + 111.051 (T/100) - 35.6543 (T/100) 2 - 

32.7 365 (T/100) 3 + 11.9649 (T/100) 4 

s f (T) = 0.712406 + 1.63522 (T/100) - 0.22175 (T/100) 2 + 

0. 0763933 (T/100) 3 + 0.02347 (T/100) 4 

s (T) = s.(T) + (h (T) - h-(T) )/ (T + 273.0) 
g f g f 

r(T) = A + B (T/T cr . t ) + c/(T/T crit ) + 


D(T/T crit> 2 + E(T/T crit )3 + 




where A = 19.66, B 


E = -2.152941, 


p(T) 


" crit 


r 

e 


-15.5499, C a -11.0722, D = 9.1141 
F = 1.3127 , T crit = 405.50 


where P ..= 113.53 
crit 

2 

0.9767 - 0.03664 (r) + 0.001338 (r) 
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APPENDIX - B 

B.l COST OP ELECTRICITY 


Year 

1972 

1976 

1977 

1973 

1979 

1931 

193 3 

Cost 
( Rs/KWh ) 

0.36 

0.43 

0.45 

0.50 

0.55 

0.60 

0.65 


B.2 FORECASTING FUTURE ELECTRICITY COSTS 

Let a variable x be defined to represent the 
number of years from the present. The cost of electricity 
( c ) then becomes a function of x. Functional relation 
between c^ and x has been found in the form; 

c 0 = a/[b + exp (-cx) ] (Bi) 

If c^ represents cost of electricity for the year 
e i 

corresponding to x = x^ , i = 1,2 , . . . ,n, the sum of 

square of errors (when eqn. (Bl) is used) is 

o n 2 

E“ = 2 [c - a/ btexp(-cx^) j (B2) 

i = l e i 

2 

The constants a/b y c are to be so chosen as to minimize E . 
This has been done using the general minimization tech- 

o 

nique, considering E as a function of three variables 
a,b and c. 

The following results has been obtained: 

a = 1.293 
b = 0.97986 
c = 0.09338 
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APPENDS - C 


Rofor corrrputer programme and data file 
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TABLE 4.1 

Variation in optimum insulation thickness with inside tempe- 
rature and sink. * 


THERMOCOLE INSULATION SUGGESTED 




DESIGN 

APPROACH 




Inside 

T emp .0 ° c 

Sink 

(m) 

Q, 

des 

kJ/nr 

Total cost 
Rs ./year 

Optimum insulation (m) 
t t R 


0 

109163.62 

31147.23 

0.0346 

0.1472 

0.0626 

10.0 

3 

95396.20 

30239.91 

0.0784 

0.13563 

0.0572 

6 

347 33.99 

30124.34 

0.0693 

0.1202 

0.0503 


9 

75050.67 

302 33.10 

0.0699 

0.1209 3 

0.0510 


0 

1198 31.3 3 

35125.94 

0.0978 

0.1712 

0.07 34 


3 

103846.9 3 

33801.41 

0.0891 

0.1550 

0.06593 

’ 4 .0 

6 

91733.02 

33358.41 

0.0792 

0.13702 

0.05782 


9 

81144.82 

33298.25 

0.0683 

0.11747 

0.0505 


0 

12477 3.75 * 

37772.16 

0.1121 

0.1939 3 0.08133 

o 

• 

o 

3 

107276.88 

36110.73 

0.10163 

0.17632 

0.07469 

6 

9 4 307.62 

35468 .04 

0.0898 

0.154905 0. 06510 


9 

8 3142.48 

35297 .89 

0.0803 

0.13312 

0.05782 


0 

129459.11 

40332.66 

0.1203 

0.2105 

0.0902 


3 

110511.77 

33 340.95 

0.1159 

0.20167 

0 .0853 

"4.0 

6 

96709.65 

37507.44 

0.1081 

0.18701 

0.7892 


9 

35121.53 

37 2 38 . 67 

0.0937 

0.16163 0.0679 
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TABLE 4.2 

Variation in optimum insulation thickness with inside temperature 
and sinks 


POLYSTYRENE FOAM INSULATION SUGGESTED 
DESIGN APPROACH 


Inside 
Temp. °c 

Sink 

(m) 

Q des 

kJ/hr 

Total cost Optimum 
Rs./Year t 

insulation (m) 


0 

126965 i52 

36463.16 

0.0644 

0.1127 

0.0483 


3 

115874.09 

34932.19 

0.0561 

0.0976 

0.0415 

10 

6 

94915.97 

34213.04 

0.0509 

0.0881 

0.0372 


- 9 

86988.19 

33570.46 

0.0453 

0.0781 

0.0328 


0 

149402.00 

442 32.20 

0.0719 

0.1258 

0.0539 


3 

134263.72 

41944.85 

0.059 3 

0.1032 

0.0439 

4 

6 

110620.19 

40839.57 

0.0591 

0.1022 

0.0431 


9- 

100540.6 

39540.67 

0 .0508 

0.0881 

0.0371 


0 

130578.4 

40619.98 

0.0751 

0.1314 

0.05632 


3 

115592.19 

40875.50 

0.0693 

0.1205 

0.0513 

0 

6 

1116 30.31 

40833.00 

0.0650 

0.1125 

0.0475 


9 

105536.36 

40844.7 3 

0.0583 

0.1009 

0.0426 


0 

135578.83 

40185.79 

0.1033 

0 . 1808 

0 .0775 


3 

109782.22 

38 2 04.79 

0.1061 

0.1846 

0.0785 

-4 

6 

96064.03 

37 386.61 

0.1009 

0.17 35 

0.0726 


9 

84574.02 

37136.15 

0.0968 

0.1675 

0 .0707 



OPTIMUM PARAMETER VALUES FOR DIFFERENT INSULANTS, ORIENTATIONS AND DESIGN APPROACHES 

(INSIDE -TEMP = 0°c) 
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TABLE 4.4 


Percentage reduction in total cost and design heat 
load over standard design approach. 


Insulation 

Orientation 

l cost 
Reduction 

% reduction 
in design 
heat load 

Thermo cole 

N-S, E-W 

35,52 

34.96 

Thermo cole 

NE-SW, NW-SE 

36.02 

34.38 

Polystyrene 

N-S, E-W 

30,64 

♦ 

32.17 

Polystyrene 

NE-SW, NW-SE 

30.69 

31.96 


thickness values, energy saving and cost increase for 20% and 
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IG 2 1 SHAPE AND ORIENTATION OF THE COLD STORAGE . 



SOLAR RADIATION IS CONSIDERED FROM 6AM TO 6PM 



TIME (hr) 

FIG. 2.2 VARIATION OF AMBIENT TEMPERAT' WITH TIME 
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THERMOCOLE INSuLATlC' 
STANDARD DESIGN CONDlTl 






FIG-. 3* 3 : 


FLOW DI ATJ-RAM FOR OISRJOR 'P'^NAIiTy FUlCHOlf 
MSIMOD 











TIME (hrs) 

FIG. 4.1 VARIATION IN AMBIENT AND SOL-AIR TEMPERATURE WITH time 



FiG. 4.2 VARIATION IN DESIGN HEAT LOAD WITH Si 




TAGE REDUCTION IN COOLS? 






COST /YEAR (Rs) 









POLYSTYRENE FOAM INSULATION 

0 - OPTIMUM SINK 

1 . insulation thickness 



RG L 12 VARIATION IN tOTAL CG$f WITH 






THERMOCOLE INSULATION 

O - OPTIMUM SINK 
t - INSULATION TH1CKNES 



SINK (m) 

FIG. 4.13 VARIATION IN TOTAL COST WITH SINK. 



t COST WITH SINK. 




SINK (m) 

FIG. 4.15 VARIATION OF TOTAL COST WITH SINK. 








HKUtiKAM LISTING 

******************************* 


MAI M D ROGRAM 
************ 

INITIALIZATION of design vector 

PRINTS OUT FINAL (OPTIMUM) DESIGN VECTOR AND TOTAL COST 


SUBROUTINE IMTPFN 
***************** 

CONVERTS THE ORIGINAL CONSTRAINED PROBLEM INTO 
UNCONSTRAINED ONE. 


SUBROUTINE POWELL 
***************** 

CONSTRAINED MINIMIZATION. 
GENERATES THE PATTERN DIRECTION. 


SUBROUTINE OFIT 
*************** 

ONE DIMENSIONAL MINIMIZATION BY QUADRATIC INTERPOLATION. 


SUBROUTINE FUN 

************** 

EVALUATES THE DIFFERENT INDIVIDUAL COSTS 
COMPUTES THE OBJECTIVE FUNCTION 
CONSTRAINTS ARE IMPOSED ON DESIGN VARIABLES. 


SUBROUTINE COPF 

*************** 

EVALUATES THE COEFFICIENT OF PERFORMANCE OF SYSTEM. 
COP IS USED IN CALCULATING POWER COST. 


FUNCTION CP-EC 

EMPIRICAL RELATIONS FOR THE PROPERTIES OF THE REFRIGERANT. 
SUBROUTINE OVAL 


*************** 


READS THE DATA FILE 

CALCULATES THE DESIGN AND ACTUAL COOLING LOAD, 


SUBROUTINE S 
***************** 

COMPUTES THE VARIOUS SOLAR ANGLES, 


SUBROUTINE HCDFTD 

***************** 


COMPUTES THE INSIDE HEAT TRANSFER COEFFICIENT., 
*********************************************************** 
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f COuPHTATIO.'iS FOR OPTIMUM PARAMETERS OF A COLD sTQRAGfe; 

■f 

}: 4. *tf*lfM!t*4''**1-4'*l* + * + * **♦**♦♦♦**»*»***♦+*♦ t*+**t»***** +; *' 


* ‘HIS PROGRAM GIVES THE OPTIMUM INSULATION THICKNESS Atr,^ Wf 

I ; i 'i !,) -,l^; Vf ' ij AN u ON THE CEILING , OPTIMUM LENGTH A5& 8PF*l— * < 
! v u n A GrvgM VOLUME. THE HEIGHT IS DEPENDANT 0» THF l.ENGT^f 

* * - yj?*2 , i55 , 2 TH * F ' 3 5 A *I^EN VOLUME.it ALSO GIVES M ... 

tuis program was developed by p.l.vemkatesh at i.i.t kan^um 

0 9/10/84 FOR THE m.TFCH THESIS. 

'•Lie oain Program 

XC l )= INSULATION THICK NE'SS (ABOVE GROUND LEVEL) 

Xl 2) SIMULATION THICKNESS (ON CEILING) 

X( 1)=SI IK 
X(4)=LE 1GTH 
X ( 5 ) = }JREADTH 

0 1 MED SI ON X (50 ) 

CO '*10!J/GRP1/R0,EPSL0N # DELTI 

C»vimoj fj 
C }0'inn/pP6/Ei 

0 AT A oPSLON, DELTA/-. 0! , 1000000.0/ 

RO=AHS(EPSLON*DELTA) 

PELT l=DELTA 

X(1 )=O.Ob I * THESE APE THE INITIAL VALUES OF THE DESIGN V 

X(2)=«),U3 

X(i)=0-1 
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X(5)=I9.0 
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1 COST 3 .COSTS, Cl ,C2,C3,C4,C5,TC) 

TYPE 608 

FORMAT (//,3X, 'THE REQUIRED CONVERGENCE HAS HERN 
i ACHIEVED* , // , 3 X , 'THE FINAL DESIGN VECTOR 

F ll ' \ I (J'fxj INSULATION ABOVE GL * \,i t 0. b/ /3X , 'f PT 

1 INSOLATION THICKNESS ON ROOF = , ^ 10 . 5 // 3 X, 'OPTIMUM sink - 5 V 
2,F5.3//3X, 'OPTIMUM LENGTH =*,FS.2//3X, 'OPTIMu" PhLADIH - , 
3F5.2// ) 


TYP'P 24'; FM 

FOR.'aTC/ , IX, 'F - MINIMUM * ',EI5.6) 

T Y pf* ^5 w X 

T 1 P li A b o U # » t C U P l C US i 1 jCOST?HcOST 3', J COST5, HEIGHT ? S I u K \ T J ft S 1 , C 1 f 
- *Aat{//{oxJ 'c6?='^O/10^,'C0ST of I NSUL AT ION /CUbJEj* * , P 8 




A m / I * A B l W 17 V I * M I JL v a s ' — # * 4 V « £ r a v * # o v * * • • * « *„ ■- ' - w r j 
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1g he a r loa5b' # | , is,I/iox, i actual heat load- J 
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1 I 
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s( r 
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C i J M 1 , ) , » / R I f l / A I J 0 S T 

C . 

C ! ’Mil) /MM, 3 /s 
A MM 4 = 0 » 0 
i:J'»>KX= t 
Ct'l i'l I 7 HO 
r ts:i 

; M 1' / 1=1, M 

Hi>( n = x c i j 
fH=r . 1 

i ■> X X = t > X 

TKQUnrx.ftq.nGO CO 7« 

\ .jPM A = A IjvJ'i 
Atjl. =ALPH A 

I fS,'!^ T 5 VcSgfl'^l:SSr^rl!e?V?c?' ME;ir ’ HT ' ST, ‘ K ' T1Kt,, ' 0Dt ' s ' t,T0 ' T ' 

Ft “F Al/t 

ALP 1 A=A)jP'1a+DX 

^ftfeT^.Sgf?^^?§;g§ffl;gS???s HFIGMT ' SI,iK " rT)iS1 . 

IFCFl.LT . frejgo Tu it 

F2=FRE 
AL2=ALPHA 
OX=3.'}*i)X 
ALP!iA= ALRHA + DX 

5SSkTSV3S*i i l?5i?8S:gS??l;gg;i25- IICECMT ' BI " ,t ' TT,, * 1 ' on “' OTOT * 

TFCFRE.LT. F23GO TO ?7 

F i-FRE 

AL3=ALPHA 

GO TO 66 

F2=FRE 

AL2»ALPHA 

GO TO 3 i 

F3»F*C 

ox*ox/ 2 .o 

AL 3 * ALPHA 
ALPHA= ALP HA-DX 

^TS?M§f? 8 ?: 8 Sta: 8 ?!??t He,BMT ' SH *- T, * s, ' ,M8 ' 0Wt ' 

T = I + t 

IFCDX.LT. O.O) GO TO dB 

SF3sfRE 

SAL3=ALPHA 

GO TO 99 

SF2f FR£ 

J? At/2 3 ALPHA 

CFCUGf. 5, XNO.DX.LT, 0.0, AND. 


.FUGO TO 111 

>5*o ■* o * 


FftE.GT.FJ ) 

“T.Ft 5 



0X 3 -DXI 
GO TO r 


IP s 


Ait2wA©*ffi4 

■■ . 

Ai.laSAU ' ‘ ‘ v ' ■’ 

KJ3.SF2 
4M=3AI 2 

4 4=1. Al3-\u1 )*(AL3-AL2)*(AL2-Al1) 

( ( . uT . 1 » o L ^2 0 ) LU Ti J )5 

H (tL2-*L1 )-(F2-Ft )*(AL3-AU ))/AA 

f'Cr.uO.'i.i )f,ij t r 4H 

M=( F3-F l 1/( AT,2-At 1) -C* { AT 2+AL1) 

A = rr-i,i* A!. t-C*ALl * At l 
‘.LS' r = -rVC 2 .A*C) 
i 4* n * 4 LL'i'+C * Al. ST* ALST 
1 i/’H As A u.i T 
! /= l t 
i 1 > 1 . X = l «HL>. » t 

I'M ( ,Cu°.CuSTl .r;./iT2, HEIGHT, SINX.TINSI .ODCS.OTOT, 
1C 1 , Gi'-ITM ,C t ,C2 ,C 3 ,F4 # CH , TC ) 

} ~ l/'.F (. ( P -v ' ) /f *-* } 



■; v • ; t / > * I i’U 

^ tt ) •G , UF2 )Gu 10 } J3 

i i j 1 - A j 4 2 


§*' i = t 2 
I* 4 = ! ii 
•i" J TO $ti 

li ;, f Fn.GT.F2) GO TO 12? 

ALi = AL2 
A 1,2 = ALST 
Pi = F 2 
F/ = P 4 
G i T ) 66 
Ai., 1 = AU<5T 
F3 = F i 
<■ J TO 66 
6;, 1= ALST 
Ft -p i 
( Tj '■’ ) 6 t> 

J£^.fi s I* L, £x $ L2><i0 T,J 13 

{f'Cb’I.GT.b}) GO TU 14 

A 6 4 s AJj2 

Mj2=ALST 

F)=F2 

F2 = FM 

GO TO 66 

r i =pft 

Ah 1 = A L, S T 
>, < t • 6 1 > 

CF(KM.CT.F?)GO TO 15 

Aul=AL2 
F l =F2 
A62=A1 jST 
F2 = F0 
•JO TU 6u 
F 1“F0 
A!j ieAjjSf 
00 TO 66 
ALP'JAsAHST 
A 6061 = ALPHA 

£T ~ E .. 

RafURh 
FsF 2 

CONTINUE 

RETURN 

FAD 


The subroutine for evaluating the function value 


SU8R0U TINE FUN (X,F, COP, COS'*! ,C0ST2, HEIGHT, SINK, TINSI ,Qr>ES,QTOT, 
iC*)5T3rCnsT5,CI ,C2,C3,C4,C5 ,fC) 
real lencith 

COMMON U 

COMMON/ PP6/F1/B QRL'E/GJ , NCNSTR 
COMMON / GRP i /RO » EPSLON , D£LT1 
COMMON /PP*/GJSUh/T/OP/$m3/S/T2/AI.M» A 
DIMENSION X(50) ,S150),GJ( 400) ,00(50) 

<>0 J 1 * 1 , N 

XCT)*()0|l|±AIiPHA*sa3 

Mm i&JtcoI{2 > 

v Ijf Tt MS ^ESs t ■'■■■'■' ' 

f* ; " n - u •' - 

; 'I'Svi ; 

- JMIir^ MJH&cryvsu.;,,-. , , 

X 5 -*> (5 i 

<aafoQO.O/(X##3C5l 

ME I GH tx'Ab ; 1 1 *S=X C I ) ! T 3 *<SB=X ( 2 ) i TI NS 1 *6 . 7* *X (I ) J SINK = >*( 3 3 

OP MG f HsX ( 4) ? BRFADTsXf 5 ) 

C Uu T’VAL(Xl ,X2 »>20,X3,X4,X5,X6, ODES, QTQT ) 


r \.M t = 'XKr 1 • rosil TS (HE INSULATION COST BS/CUBMET 
r ♦.»n = 5 « , 6 h,o n COST 3 IS COST OF RfcFkIGEP.ftTI. 4 G MACHINERY * 8 /KM 
C.»r,T5=3«». •• COSTS 1 ? ' *' " ■ * < •' / ' C\ 6‘ 

f .75 it h? is t> ■ ■ * ■- r r *. fT 121 

•■'•-t .2 11 SF IS fhfc DESIGN SAFETY FACTOR 

P-U. 1 ’ 11 h IS THE RATE uF INTEREST 

L= 2 C 1 I l, IS the mfs OF the project 

Ai.=f. 

’ .!-K i . )MH .C)/(P*f l.C- + R)**tL-l ,0)) ! IPNF*PRESC«T WORTH 
r!sCOGTl*C?.n*?rs*(HETGHT-SI :)» T . •’ ,* ~ r • ’ r " t *'\ * $' •>'" '•! * 

HU M G T i i + s .» k k A u T ) + L E <GTR 4 BREAD T *TI HSR 

/. Hi :oH Lt/*(.TM * h» i. A0T5/AT, 



I = u 


+ l. } *■ K Kl '• 


1 ’ • = c o S T 5 * b L: (V G T H * B k F A D X * S I N K * i . 2 / A L 

********* ***♦♦* r-1**4 ♦♦*****.*'*♦***♦* #***4:*4t*t ******** 

i + F2+C3+C4+C5 i * FI IS THE OBJECTIVE-TOTAL CoS f 


.1 INSULA TT UN 
II INSULATION 

I i INSULATION 

I I INSULATION 
I LB) 


THICK wESS 
T ; 7. C : < -TSS 
THICKNESS 
THICK WESS 


ABOvf CL(UR j 
ABOVE (~L(L!U 
ON HOOF (J., n ) 
OH HOOF (UP) 


C - E' 1 

} k * ;< * : 

) C I ) - * * » 0 H / X f 1 1-1.0 
c.icnsxct )/o.is-i ,o 
U.J( 1 ) =0 .0 7/X ( 2 ) -1 .0 
U K4) = x (2 J/O, 20-1.0 
.) (s )=-xni ; i sink 
C.iU,) = X(3)-8.0 11 SINK ?UB) 

C?) = l b,u/x( 1 1 LENGTH (LB) 

• r (B)=XC4j/30.u-l.O I* LENGTH (UB) 

' ('/)~1 ! J.0/X(5)-1,0 II BREADTH (LB) 

u0a«,') = X(5)/30.0-i,0 II BREADTH (UB) 

* * 1* ****************** *** ****** ****** ****** ******** ******** 

iC iSTH=10 
GJStiMsO.O 
f' t 3 1 f 1 = 1 ,.»CNSTR 
TF(uJU) .Gr.EPSLPrOGO TO 30B 
U.ESU AsGJSUM+1 ,0/GJ(I) 

•Ti TO 310 

G,lSn.M = GJSUM+C2.0*P.PSLUN-GJ(I))/EPSLON**2 

Cu.i 1’INUE 
F=FI-RO*GJSUM 
RE TUB J 
E. 4 0 

* *.**,* 4 . ************:******* ******************************* 

S a B R 0 u TINE C 0 P F c CO P ) 

I;*************************************#**************** 

TC s 4 2 . 0 ?TE=-8.0;TUC=5.0; TSH*5. 0 f EFFC*0 . 8 5 

I S Ta R f ”^LCU LA TE S^fu ?”(!()?"? n fT CONVERTING KW-Hfes' TO ToU-HRS 

•tT=te~ " ‘ " 

T5=TE 

TlHrTEETSH 
T i=TC 
T il;=TC 
(4= I'ID-TUC 

PaPB(TC + 27.3.15)/PRCTEt273.15) 

H1='IG(T1 ) 

iUO*!U fCPCTl )*(T1D-I1) 

Sinlsi fCE(Tl )*Af J OG((TJD+?73.0)/(Tl*273.0)) 

S20sS10 
*>3»5G( T3 ) 

T2D*(EXP((S7n-Si)/CP(T3)))t(T3+273.0) 

T2HBT2D-273.0 

H2f)«MG(T3 )+CP(T3) * (T2D-T3) 

H2»(H20-H1I))/EC(P)+H1C 
H4D»HF (T4l>) 

H4-HFC 1 4 ) 

H5=H4 

COP* ( H1-H5 ) / ( H2-H ID) 

RETURN 

E! # D 

********************************************************** 

**************** 

%mm mu«$ THE EFFECTIVE ELECTRICITY (POkkFH) COST, 


liSHtTS"’' » » . • 

5:1 5a " 

ZontUU ' 

r ** rt 


, 


CvgllMlfc 

al*L 

CEFT Vh*( ,U 
uU 3u0 1=1 


tt^RB+EXPC-O. 09338*1) ) 


: 


iU 3u0 1=1, ( 

7K*'l Vh=CEf T¥F*CL( 13/(1 .0*H) **( l-l ) 

' )U7 1 *.UE 


' i ST 2 =C E F T * t / A 0 

t.. T 1 J !J f 
S- I; 

M ************* ********************************* 

; "t-T-ICAL Fii.'iC l'IO\R FOR PROPFHTIES OF AMMONIA 

* 4 V, 4, ** % t4»tt*»*»*44*4***»****** **************** 

* <: n rx, c : ( t ; 

T F ( 1 . i ' ,f . 1 • * . 0 1 C.n 10 

• = , 341 .fM 4 r-i 3* (T -1^.0 *0.050 15 * (T- 1 (‘ « 0 3 * *2 ) 


= 2.1 1 ’ 7 5 + 0 


Vi b qiH*( i*44S.C*0.UO('C6OCb3b*(T+4r»,()**2) 




! l» ! '( l»F ( T ) 

/i oi, .0 

i;k.ftU5 + 462.25*A + 28.7J6B*A*A + U.836*A**3+5.:*621«+A**4 


/lvu .0 

1 1 *3.36+1! Ue51*A-&5.6543*A**2-32.73b5*A**j+lJ .Of,49*A 


&r (T) 

06 + t ,6d52? LA-0.221 75*A**2 + Q.07b3988*A**3 + C» .02346 


'■••fiiiCTlOu 6G(T) 

•H = ■ S F ( I ) + ( 1 i G ( T ) - H m ) 3 / ( T + 2 7 3 . 0 ) 


F.HCTIOSV PR ( T) 

roi r=n,s.so 
PCRIT=I1 3.53 
Asl’i».b6 
B = -i5. 5 4-0930 
C=-l 1.07219 
I>*9.llil07 
G=-2. 15^941 
F=l Jl2o9 

itF*(U t S!-l/TMlT}**l / 5?/?T/Tr!9? /TCRIT)f * 2HF#UT/TCW,T) ’ f3} 

PR=pcrahsxp(rtL) 

RETURN 
R (!) 

FUNCTION KC(R) 

RCsO.. 9T6695*O > 0366.432«Rf0,0Ci 337988* P**2 

sud 

♦..Me*****-*************************.**********#************** 

*************«********************+********♦******»*»******«*** 

SUBROUTINE y VAL (T INS r TXNSK, TINS! , SINK# L tHQTh ,BPEAUT # HEIGHT f yD£S 

♦ &** ******* * ********* *******t**************,i« ****♦*#****#,#♦♦**** 

IN TEGEH Oft I ENT , TIME ( 2 4 ) 

PEAL. IRAD(24),IMB(24) / iNW(24),INSt24) e INN(24),IMP(24) 

DIMENSION TEMP ( 24 ) .SOLE C 24 ) , SOL* ( 2 4 ) ,$0bN < 24 ) , SOLS (24), SOLE ( 24 ) 
DIMENSION T(6),0I«F(24),T0INF(24) 

Hi- \h KINS# K. 1 6 ) t LENGTH , M AXIMP,NCH , IDIFFC24) 


OPE ■ C NIT S 2 2 7FlLF“7jr:DAT T T 


_ _ ... WD»* 1 / « i J 

—INPTJT STAPfS-— 

ALAT, MONTH 
TlfHTxPHt 
NLA YEP 

(T(U,K(IMtJ .NbAYER) 

FINS 

TE'IP(I) ,iPAD(I),IDIFF(I) 


READ (22,* 
AS AD (22,* 
RE AD (22,* 


a;,*** -L, 7 * 3,1 417/itTM 
URANGl,s*90.0 
DO 5000 1*3,24 
O’-.Ai ■ “OP AH j U 

i! ;i*,Ds.lP4 •■Gl ) * J. 1 4lWh/i 8<) 0 

H ■» \ 1 )=l 

VC I i. ( n.GE.b.A,vr.Tl»‘FfI).LE.l8) GO TO 25 

; pi i )=< .< 
r « • T 

! ■■ { f is: ... 

! *> t T 1 -r < . i 

T '■ '•>( 1 J=(' 

;•*. to 175 

■"••< SOI, A.’P 1 .HVAbjKDuTH, BETA, GAMA) 

V" • "" D 3 S C *, A ‘ ) 

j.G-. .1 1.1415926/2, U))G0 10 1? 



/ 

I. Vi ? x ) 

V.. 

1 

w*' f *« 1 

~0^ 1 

w*»h 1 

.1 » 

4 1m fc‘* ^ t S !K 

V* 

roc Ai*f 
C OSTHk 

00 to 

1 7 

c 0 ST M f* , 

1 i 

r*v:C(IJ' 

c 

r i«' ( t r :*i 

Ab^ i A Av 
G 1 TO 

<1 h 

A IjPb \'d< 

4 7 

TP( AbP' 
r to 1 Hr: 
GO 10 

,1 *> 

COSTCO 

4 i 

T : OJ): 


I ■' U ) = U H Ai> ( X ) + TOIFFC I ) ) *COSTHS 


K / 
A , 


5 7 
s 

c 


p i 
173 

175 


c 

c 

c 

c 

c 

c 

c 


315 

316 

If -I , * ■ . 

31« 


i r 


i i<i 


i .1 : > 


= Aas C G kn A - 3 . I 4 1 59 26 /2 , 0 } 
3o 


^ ;* * j • v / • * ^ i j 7 & f & » v 

:COo(BFTA)*COS(AT,pHAF 

:0 * 0 


)G0 TO 3? 


4 7 


(D.GT.U)GO TO 46 
A Bo C GAM A -3 . 1 415926) 


ABS(f»AMA + 3, 1415926) 

AG . GE . ( 3 . 1 4 1 5926/2 . 0 ) 3 GO TO 42 


:0 . (; 


WsABS(GAMA+3. 141 5926/2.0) 
EriAW.GE, (3. 14 15926/2. 0))G0 TO 57 
W=C05C BETA) *COS( ALPHAW) 


T F ( 1 1 HE ( I ) • LT . 1 2 ) GO TO 52 
AGO.iAKaAB5Ci,ArtA-3. ! 415926/2.0} 

an to 56 

A ijP<i = ARSCGAtoA + 3. 1 41 5926/2,0 

ft-' (ALPHAS, 

COS’lfiW _ 

GO 10 55 
CJSTH#=O.U 

r) = (TKAr>( I) + lDIFFn))*COSTHW 

l!^(T) = CIPADC n + n>rt-F(T))*SlN(BETA) 

if It ime c iIIlt.!/ Igo^tcGsI'" 

HRA;!Gl.=HRANGL+15.0 
Gn TO 171 

6HANGL=rtRANGL-15,0 
RKT A 0= BETA* 1 80.9/3. 1 415926 
C.ViAD*GAMA*18Q. 0/3. 1415926 

l )-TFMPC l ) t ABSPTY*IKEf I )/Hf) 


s:jl,s ( 1 )=TG«iP f 1 ) + ABSPTX *1 NS (I)/HU 
SOLw ( 1 ) »TEMPC X ) + ABSPTlf *f N* (I) /HO 
TEMP ( X ) 4 ABSPTY )/H0 

TEMP ( I ) + A DSPTY * 1 WK f t ) /HOP 

clTIZbJlWjroF^^h^TP^~J"I~DS'. 


salmi) 

Snop(i) 


C 0 H T IhOE 

Sul“”fTON“ 5 f"BtE 14 T ANDCSOLAt^ ' £ 4 r Afct$ 


HERE THE TEMPS ARE SUMMED UP FROM 10 HOURS m 15 HOURS, 

SlGflSo^)' — ' “* " — ^ — - 

no 315 1*10,15 

SIGTA*SIGTA+(TErtPU>-TT) 

CONTINUE 
S|GTS*0,0 



70 319 1*10,15 

s i rr i’=si G?tu c sour t i ) - terp id) 

C }•* I r .ml 

r -i T 6 T A = i .(• 

r Vi l u r . I = » ,24 

I* HOT =r,SlGl a + ( lFRP(I)-TT) 

r i i r on 

fol , 1 • > = »> 

! i 1 1 o i _ i , i \ 

1AI-, 1 r.s «•> 16 rsi (SCLS(t)-TK-iP(I)) 
!V I ( H. 
r .i ■. 1 1 =« .< 

' • 1 U7 T = 1 ,2\ 

r,.f - i s tori- ■♦(.s'MjF tn-ir^u )) 

r t -T I n-: 

IAS , r 



1 > i 


1 J ' : 

c 


1 V 1 


i 1 :t '/ 

C 


9*7 
•illl 2 

c 


ts»37 

c . 

1 J 11 


c 

c • 


r > r ' j t , i ™ t / , ? 

■*•*' I 5 19 1=1 , ?4 

rs ! n = ISIGTH+ t SOLP c I J-TEMPU ) ) 

C • « >' r l 

— .__SUM*ATIGP FNDS. 

7 A rToy "”u T H tR MAlTi^S fSf AF err' 

r • i h .-5 j = o . u 
,v ’ ’ <s i 1 = 1 , ,LAUP 
IMS I ='i’HPST+t lf)/K(l) 

C-» *1 i.-.im*: 

■1 AX r n’=TE..P( I h ) > 

M< r i v p ~ 9 A >' t l’E.'.P ( 1 4 ) , 7 EP P ft 5 ) , TEMP (1 6 ) j 

'Hi Mi 

M = $ R E 4 [) T 
M = H£IGIiT 
S=S l i K 

CA.GC.'f.A.rrnri of ipfiltepation load. 

r ? J ™ i • ^ 0 ^ j 
HMFT=0,C 

7=A*B* (H-S) 

1 ) + o * r 4 J708* ( ALOGC V ) ♦M ) + ^ 1 . 4 2 6 ♦ C ALOGt V 32*( ALOGC V ) **3 

Dn 'ill 1 = 10,15 } } 

A A= l . J 526-5 
00=4 , 7076-9 
C = T K -1PC i J+27 3, 15 
0=6-17, 51/C 
! ■: = C C - <4 1 C . CM ♦ * 2 

PJ=22J .22 B/EXPf (7.21379+f AA-BH*C)*n*CP-l .0)} 

VA=2 9 27 . 9066/ ( (1 0132-PVm .065) 

5 of Sr * 1 1 ?©f 2rT C *** ***■ 

CU'-TT:10E 
T v >IMFT=0.0 

VsA*t)*(4-$) 

|CH = i J.O^C8-54,3906*A|UGCV) + ll ,426*(AL0GCV)**2)-1.132*(ALPG(V)**3 

00 iVl? 1=1,24 0 V 
AA=1 .152E-5 
DB=4 .787E-9 
C=T6MP{ I)+273.t5 
0=647,31/0 
6=(C-21O.0)**2 

PS=221.228/EXP{(7.2l379+(AA-R0*C)*E)*(D-i.O)) 

PY = t?Hl*PS 

VA = 29,27*C*9 J 80bfc/((l,ni32-iPV)41,OF5) 

r.jr /F(I)=A*rt*(H-S)»ft!CH*CP4(TEMP(f )-f J )/( VA*24 ,0 J 

Tin^Fr=TQINF(I)+Tt>nFT 

CONTINUE 

CALL HCOFTDITI , MAXTMP , THPST ,T1NJ>*K INS « Hf),U) 

CALL HCOFTO ( TI # ttAATMP. THPST .TINS*. KINS, HO » UR) 

Ut=l ,0/Cl .0/2, 0tl,0/H0+TWRSi+TINSl/KTNS3 

i A5t0*A*(H-S)*(SIGTA»?.0tSIGTS)+U*B*(H-S)*C2.0*SIGlA4SIGTE+SlC.TH) 
♦UR*A*H»CS1GTA+SICTR) 

2= l*2.i.*F>*(01-Tl)*(A+S3*5.0 

O-S >. ( 2 .C.TSJCT-A.TSIGTr* 

fex|i^m T ?? n * ,A ‘ fl, * 24 -° 

tfu|«N “.'M , . 

* 'T / ^ * 1 ^ 

T? , : .„ T - r: ._ w „ rT „ rTrr , 




■ who 


p IV 

If < 

i f r 

f f. f 


1 ■ 

fa J )» 


! i- 
! r 
1 K 

i i- 
lb 
It 
I b 


At 

*ZQ* 

• H 0 , 

Q, 

.KO, 

• Kf>, 

,Kv. 

r ^ 


i ,f. 


,i.i 
: < * ■ V 


. E 0 , 
.M, 
.LO, 

= CP- 

? v r 




1 ) PEC=-21.483 

2) PEC*-I 3.7SP 

3) '»tC = -i. 7675 
i)i cCs7.P’.it> 

) ^-C = 17 , i«l / c 
0) 06C=2 / . 0? 

7 )t'cC = ?1.77 
b ) t *’,(' = 1 b , i * 

0 )Sm : .C=M 4 375 
At jOCC =-6,6325 
1 1 ).:- C = -1 ( .( 3 
J 2 )l^C=-?/,65 
,14 1 7 /M .0 

>n S t C ,^0 ( CkAJj ) ) 4COS C DFC ) 



r 

c: 


1 ip i 


‘ n o X V o i v b v i 


C 

? i; t -• * > 


t 1 


* C‘ 


j .u !. ; . m rih hcof th ( s tt , m axt , sthpst , stink , kin , sun, si ! ) 

U ,v A T , fj| it XT f F 1 N 
i ’ ,/i i=c.< i 
1 l =5 * 1 + Dfc.tiTl 
• i f - 5 . IN* (DEI.T1 

s !=! .o/(l ,0/HL+l .U/SHU+STHPSTtSTINS/KIN) 

>' i >M=So*(.‘HXT-sTl) 
ni LT2 = FLUXt /HI 

•■’."•-U = ABS (. PELTl -DFJLT 2 ) * 1 00 , 0 /OhLTl 
T •’ ( pi*;KIi 1 F » T/r . i . f : )G 0 TU 1000 
Of J,T 1 =QELiT 3 + 0.00025 
i>i'TiTl=DELT 2 
Gii 'Xu 70 00 
R >•, T OR! ; 

♦Mt*******^********************** ************************ 





f 

* t # 

7 » \ / , * ,,-.22,0 

1 f * ’ 

b i . 



4 f 

i . , i, , I .0 

p 1 ® 

• 9 

o . o , i - , h 

*' / '» # 

4 9 

1 . 1 • , 4 

*' 1 « 

1 # 

* '■ 9 4 , f ' W 1 

i, s 

/ , 

• ^ , * ' . ^ 1 

#1 u • 

; 

1 V J l . 5 1 , B 2 , 6 0 

■/'■ ' 4 # 


'• 5b. i22 

J 1 *-■ 

ti ? 

v,7M. . 1, JtfS.lfl* 

1 j . 

-3 | 

2'»fs>J.79,2ii( ,2b 

■* K 

? ...» <B 

r 

Si f;9. h2 ,20#, 37 
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